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摘 要 在 心理 学 、 教育 学 和 临床 
心理 、 临 床 效果 等 随时 间 而 产生 


在 短 时 间 内 对 个 体 i 


旺 等 的 动态 变化 及 其 作用 机 制 。 近 年 来 ， 密 
追踪 的 研究 分 析 仍 停留 在 较为 传统 的 方法 。 方 法 学 领域 已 涌现 出 较 多 | 


中 山大 学 心理 学 系 


医学 等 领域 , 越 来 越 多 的 研究 者 天 


aT F 始 关注 个 体内 部 的 行为 、 
的 动态 变化 ,重视 针对 个 体 的 差异 化 建 模 。 密 集 追 踪 是 一 种 
行 多 个 时 间 节 点 密集 追踪 测量 的 方法 , 更 适合 用 于 研究 个 体内 部 心理 过 


追踪 成 为 心理 学 研究 的 一 大 热点 , 但 许多 密集 


j 于 密集 追踪 数据 分 析 


的 模型 方法 ， 较 为 主流 的 模型 包括 以 动态 结构 方程 模型 (Dynamic Structural Equation Model, 


DSEM) 为 代表 的 自 上 而 下 的 建 模 方法 ， 以 及 以 组 迭代 多 模型 估计 (Group Iterative Multiple 
Model Estimation, GIMME) 为 代表 的 自 下 而 上 的 建 模 方 法 。 二 者 均 可 以 方便 地 对 密集 追踪 数 


据 中 的 自 回 归 及 交叉 澡 后 效应 进行 建 模 。 


关键 词 ”密集 追踪 ， 时 间 序 列 ， 动 态 结构 方程 模型 ， 组 迭代 多 模型 估计 
1 引言 
在 干预 治疗 或 自然 发 生 的 事件 中 , 个 体 的 想法 


成 不 变 的 , 而 是 会 随时 间 的 推移 发 9 
究 方法 能 够 在 一 定 程度 上 对 人 群 中 存在 的 现象 进行 
内 部 心理 状态 的 动态 变化 ， 也 不 能 详细 解释 事件 对 个 体 的 情绪 、 行 为 等 产生 的 持续 性 影响 
能 导致 研究 者 对 干预 效果 及 其 潜在 机 制 产 49 
Setodjietal., 2019)。 因 此 ， 为 了 更 深入 地 从 个 体内 的 角度 研究 变量 间 的 关系 ， 研 究 者 往 
密集 追踪 方法 , 更 密集 地 测量 个 体 状 态 发 展 和 变化 的 纪 


i 


(Setodji et al., 2019), 可 


、 情 绪 、 行 为 以 及 生理 功能 等 往往 不 是 一 


FE 动态 变化 (Vallacher et al., 2002)。 传 统 的 追踪 研究 或 横断 


密集 追踪 数据 有 测量 间隔 短 、 频 次 高 的 特点 , 常常 被 ) 


述 或 解释 , 但 难以 揭示 短 时 间 内 个 体 


E 错 误解 读 ( 张 银 普 等 ,2016; 


往 需 


节 (Trull & Ebner-Priemer, 2013). 


j 来 对 个 体 的 行为 以 及 状态 随时 间 


变化 的 动态 过 程 进 行 更 精细 的 测量 ， 测 量 内 容 包括 情结 、 人 格 、 认 知 和 行为 模式 等 方面 


(Hofmans et al., 2019; Lazarevic et al., 2020; Trull & Ebner-Priemer, 2014)。 根 据 Bringmann 等 


人 (2017) 的 定义 ， 密 集 追 踪 数 据 是 指 研究 者 在 较 短 的 时 间 
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习 家 上 自然 科学 者 


点 的 变量 信息 。 通 过 对 变量 进行 重复 测量 , 研究 者 可 以 直接 获取 有 关 被 试 内 的 变量 信息 ， 实 
时 地 记录 被 试 随时 间 变 化 的 动态 过 程 (Bolger & Laurenceau, 2013)， 进 一 步 解 释 个 体内 部 心理 
状态 和 行为 的 动态 变化 等 ( 张 银 普 等 , 2016)。 

近年 来 ， 随 着 手机 等 移动 通讯 设备 的 发 展 和 使 用 ， 研究 者 可 以 更 加 方便 、 频 繁 地 收集 有 
关 被 试行 为 、 所 处 环境 、 心 理 和 生理 状态 等 方面 的 数据 ; 数据 收集 过 程 的 便捷 化 以 及 收集 场 
景 的 多 样 化 等 , 也 都 有 利于 研究 者 探究 关于 不 同 场景 和 不 同时 间 点 的 问题 。 如今 密集 追踪 数 


据 已 经 被 广泛 应 用 到 心理 学 的 多 个 领域 (Chen & Zhang, 2020; Windt et al., 2018), W, RIIT 


HY 


为 (Lanza et al., 2014; Shiffman, 2009; Weinstock et al., 2017)、 压 力 症 状 (Chun, 2016; van De 


Schoot et al., 2015)、 心 理 资 本 (Malmberg & Martin, 2019; Sturgeon et al., 2014)、 饮 食 障 碍 问题 


(Munsch et al., 2009)、 亲 密 关 系 (Howland & Rafaeli, 2010)、 人 格 障碍 (Trull et al., 2008; Wright 


et al., 2015, 2016) 和 语言 使 用 (Lazarevic et al., 2020) 等 。 


Trull 和 Ebner-Priemer(2013) 指 出 , 密集 追踪 方法 解决 了 很 多 回顾 性 调查 和 实验 室 方法 的 
局 限 ， 如 潜在 的 记忆 偏差 、 对 个 体 日 常生 活 场 景 的 还 原 度 不 足 等 (Bolger & Laurenceau, 2013; 
Sened et al., 2018)。 考 虑 到 密集 追踪 方法 存在 的 诸多 优势 ， 国 内 一 些 研究 者 开始 使 用 密集 追 
踪 数 据 来 验证 实验 假设 , 如 探究 干预 方法 对 个 体 心理 和 行为 变化 的 影响 ( 陈 明 瑞 ， 周 萍 ,2017; 
郑 文 倩 等 ,2019)、 环 境 因素 对 个 体 身 心 状 态 的 影响 ( 安 媛 媛 , 徐 慰 ,2017)、 个 体 在 日 常生 活 中 
的 压力 应 对 ( 邵 华 等 , 2019) 以 及 个 体 在 休闲 活动 中 主观 幸福 感 的 变化 等 ( 任 术 等 , 2010; 张 峰 ， 
耿 晓 伟 , 2015)。 这 些 研究 多 采用 较为 传统 的 多 层 线 性 模型 (Multilevel Model, MLM) 对 密集 追 
BARRETT (Be, REL, 2017; 邵 华 等 , 2019; 张 银 普 等 , 2016)。 然 而 ， 密 集 追 踪 方 法 
的 优势 更 多 地 在 于 探究 变量 在 变化 过 程 中 的 动态 效应 ( 即 变量 随时 间 的 变化 可 能 受到 变量 本 
身 或 其 他 外 部 因素 的 影响 ; Hamaker & Wichers, 2017; Molenaar, 2004; Molenaar & Campbell, 
2009)， 因 此 对 密集 追踪 数据 的 分 析 方 法 提出 了 更 高 的 要 求 。 

为 了 帮助 研究 者 们 增进 对 于 密集 追踪 数据 方法 的 了 解 , 下 文 将 首先 对 常见 的 密集 追踪 研 
究 方法 及 其 优势 进行 整体 梳理 , 然后 介绍 密集 追踪 数据 中 较为 常用 的 分 析 方 法 (MLM), 以 及 


目前 较 新 的 分 析 方 法 , 主要 包括 动态 结构 方程 模型 (Dynamic Structural Equation Model, DSEM) 


mn 


和 组 迭代 多 模型 估计 (Group Iterative Multiple Model Estimation, GIMME)。 为 了 更 直观 地 比较 


MLM, DSEM 和 GIMME 在 实际 建 模 过 程 中 的 特点 ， 本 文 分 别 利 用 这 三 种 模型 方法 进行 实 
E 的 演示 和 比较 。 最 后 ,为 了 帮助 研究 者 针对 自身 研究 的 实际 情况 灵活 选择 适合 的 模型 方法 ， 
本 文 还 将 对 各 种 主流 模型 方法 的 优 缺 点 和 适用 范围 进行 总 结 ， 以 供 参 考 。 
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2 密集 追踪 研究 方法 及 优势 
2.1 密集 追踪 研究 方法 


密集 追踪 方法 通常 又 称 动态 评估 法 (Ambulatory Assessment, AA; Bolger et al., 2003; Sened 


et al., 2018)， 其 另外 几 种 常见 的 术语 有 : (1) 经 验 取 样 法 (Experience Sampling Method, ESM; 


Csikszentmihalyi & Larson, 2014), 广泛 应 用 于 人 格 和 社会 心理 学 领域 ， 主 要 使 用 主观 问卷 进 


行 测量 ，(2) 日 记 法 (Daily Diary; Bolger et al., 2013) 和 日 常 过 程 研究 法 (Daily Process Approach; 

Tennen etal., 2000)， 通 过 被 试 的 每 日 记录 ,追踪 有 关 变 量 随 时 间 变 化 的 过 程 ; (3) 生 态 瞬 时 评 

估 法 (Ecological Momentary Assessment, EMA; Stone & Shiffman, 1994)， 常 用 于 临床 和 药物 研 

究 。 典 型 的 EMA 研究 收集 被 试 在 日 常生 活 中 每 日 特定 时 刻 的 状态 和 行为 数据 ， 在 这 些 特定 

~ 的 应 答 时 刻 中 ， 被 试 或 是 被 动 接受 来 自 外 界 的 提示 ， 或 是 自主 发 起 某 些 行为 (如 ， 被 试 需要 
在 感受 到 对 成 六 药 物 的 强烈 欲望 后 进行 作答 ; Lazarevic et al, 2020)，(4) 动 态 监测 法 
(Ambulatory Monitoring; e.g., O’Brien et al., 2000)， 被 试 需要 佩戴 特定 的 便携 式 仪器 ， 对 特定 
的 检测 指标 (如 
一 般 来 说 ， 密 集 追 踪 研究 方法 可 以 被 广义 地 分 为 三 类 (Fisher & March, 2012; Trull & 
Ebner-Priemer, 2013)， 包 括 自我 报告 式 的 自我 评估 ; 观察 式 的 自我 评估 生理、 心理 或 行为 
ms 的 自我 评估 。 在 自我 报告 式 的 评估 中 ， 根 据 不 同 的 实验 目的 ， 被 试 可 以 在 固定 /随机 间隔 的 
二 时 间 点 进行 自我 评估 ， 或 者 在 特定 事件 发 生 后 进行 作答 (Reis et al, 2014) REW ARIRE 
© 方式 各 有 优势 和 不 足 : (D) 固 定时 间 间 隔 的 自我 报告 方式 要 求 被 试 在 约定 好 的 若干 时 间 点 回 
答 相 应 的 问题 (Fisher & March, 2012)， 其 数据 收集 的 过 程 不 会 突然 打扰 到 被 试 ， 数 据 收集 方 
式 本 身 简单 易 操作 。 然 而, 固定 时 间 间 隔 的 自我 报告 也 可 能 因为 只 考虑 固定 时 间 点 发 生 的 事 
件 ， 而 忽略 了 其 他 时 间 点 的 事件 可 能 产生 的 影响 ， 同 时 ， 当 测量 时 间 间 隔 过 大 ， 也 可 能 使 得 


Seanet 


:) 进 行 密集 测量 。 


自我 报告 中 ， 被 试 需要 不 定期 地 配合 研究 者 完成 数据 的 收集 (Fisher & March, 2012)， 因 此 可 
以 随机 抽 到 一 些 具 有 代表 性 的 经 历 ; 但 是 对 于 被 试 而 言 ， 每 次 数据 收集 的 时 间 是 未 知 的 ， 这 
种 不 确定 性 可 能 对 被 试 的 日 常生 活 带 来 较 大 的 影响 和 人 负担 (Fisher & March, 2012; Hofmans et 
al.,2019)。(3) 关 注 特定 事件 的 数据 收集 方式 ,是 指 被 试 在 发 生 特定 事件 后 主动 配合 研究 者 完 
成 数据 的 收集 (Fisher & March, 2012)， 其 优势 在 于 可 以 收集 到 与 研究 者 感 兴趣 的 事件 相关 的 
数据 。 然 而 ， 这 种 方法 收集 到 的 数据 取决 于 被 试 对 事件 的 反应 ， 因 此 研究 者 需要 事先 对 被 试 
进行 训练 ， 使 其 能 意识 到 特定 事件 的 发 生 并 主动 及 时 向 研究 者 报告 (Fisher & March, 2012; 


Hofmans et al., 2019). 

另外 ,与 自我 报告 的 评估 不 同 ,观察 性 自我 评估 和 生物 、 心 理 或 行为 的 自我 评估 主要 收 
集 客观 数据 和 指标 (如 ， 被 试 所 处 环境 中 的 声音 、 被 试 运动 的 速度 ， 以 及 心率 等 生理 方面 的 
指标 数据 )， 而 不 是 主观 报告 的 数据 (Trull & Ebner-Priemer, 2013)， 因 此 很 少 使 被 试 产生 心理 
负担 , 同时 其 数据 收集 也 不 易 受 到 被 试 的 主观 影响 。 具体 关于 不 同 密集 追踪 方法 的 使 用 可 以 
参考 Fisher 和 March(2012)、van Roekel 等 人 (2019) 的 文章 ， 以 及 本 文 关 于 密集 追踪 研究 中 的 


实际 研究 示例 。 


2.2 密集 追踪 研究 方法 的 优势 


相 比 于 传统 横断 研究 的 方法 , 密集 追踪 研究 方法 的 优势 主要 体现 在 : (1) 与 横断 研究 中 国 
顾 式 自 我 报告 的 方式 相 比 ， 密 集 追 踪 方法 对 变量 的 测量 更 及 时 、 更 接近 被 试 的 真实 经 历 ， 可 
以 更 好 地 减少 测量 过 程 中 出 现 的 记忆 偏差 (Stone et al., 1999)。 因 此 ， 所 产生 的 密集 追踪 数据 
对 于 事件 的 描述 往往 更 加 准确 ， 更 能 代表 被 试 的 真实 反应 (Sened et al., 2018)。(2) 与 横断 研究 
中 的 实验 室 方法 相 比 , 密集 追踪 方法 提供 了 多 个 时 间 点 和 场景 的 测量 , 使 得 研究 者 可 以 探究 
个 体内 部 心理 构 念 随时 间 的 动态 发 展 及 相互 作用 ( 张 银 普 等 , 2016; Trull & Ebner-Priemer, 
2014)。(3) 密 集 追 踪 方法 的 测量 随 着 被 试 的 日 常生 活 情境 进行 ， 不 仅 局 限于 实验 室 环 境 ， 极 
大 降低 了 社会 期 望 以 及 被 试 自我 暴露 意愿 对 实验 数据 的 影响 , 因此 有 更 好 的 生态 效 度 和 外 部 
效 度 (Sened etal., 2018)。(4) 密 集 奶 踪 方 法 对 多 变量 重复 测量 ,允许 研究 者 分 别 探究 个 体 间 和 
个 体内 的 变异 (Curran & Bauer, 2011; Sened et al., 2018)。 许 多 心理 学 研究 的 理论 和 假设 主要 
针对 个 体内 部 的 过 程 (Intraindividual Processes)， 但 其 数据 分 析 的 本 质 往 往 是 对 个 体 间 关系 
bond (Interindividual Processes) 的 探讨 (如 ，Lau 等 人 (2016) 针 对 个 体 间 自 变量 水 平 变化 对 因 变 
影响 ， 得 出 “孕妇 对 压力 的 积极 感知 会 促进 其 积极 应 对 的 策略 ， 进 而 提高 身心 健康 ”这 一 个 
体内 的 结论 )， 因 此 这 样 得 到 的 关于 个 体内 关系 的 结论 可 能 是 存在 问题 的 (Curran & Bauer, 

2011). 
相 比 于 传统 的 追踪 方法 ,密集 追踪 也 具有 其 独特 的 优势 。 传统 追踪 方法 收集 的 数据 重复 
观察 的 数量 较 少 (通常 少 于 10 次 ; Jahng, 2008)， 并 且 测 量 时 间 间 隔 较 大 ， 因 此 主要 用 来 探究 
在 相对 宏观 的 时 间 尺 度 上 ， 因 变量 总 体 均 值 随时 间 的 变化 (McNeish & Hamaker, 2019; 
Nesselroade, 1991; Setodji et al., 2019)， 但 不 适合 用 于 研究 个 体内 变量 的 动态 变化 过 程 (Jahng， 
2008)。 而 密集 追踪 方法 对 变量 的 测量 间隔 和 持续 时 间 更 短 、 频 次 更 多 (Hamaker & Wichers, 
2017; Walls & Schafer, 2006)， 可 以 更 深入 地 反映 变量 随时 间 变 化 的 复杂 模式 (Bolger et al., 


la 


2013)， 允 许 研究 者 对 被 试 的 行为 或 心理 社会 过 程 随 时 间 产 生 的 细微 变化 ， 以 及 相关 的 共 变 


机 制 等 进行 研究 (Collins, 2006)。 相 比 于 传统 的 追踪 方法 , 密集 追踪 方法 在 被 试 内 水 平 的 数据 
ARZ, 有 足够 的 检验 力 检验 被 试 内 水 平 的 研究 假设 , 改善 了 传统 追踪 方法 中 被 试 内 水 平 的 


数据 点 较 少 而 统计 检验 力 不 足 的 问题 (Fraley & Hudson, 2014; Hofmans et al., 2019). 


总 的 来 说 ， 密 集 退 踊 研究 方法 能 更 好 地 保证 数 


据 的 有 效 性 和 真实 性 (Sened et al., 2018), 


允许 研究 者 从 个 体 的 角度 探究 变量 的 变化 关系 , 深入 分 析 个 体内 部 的 心理 变化 过 程 (Baltes & 


| 中 


Am 


Nesselroade, 1979; Molenaar, 2004)。 利 用 密 


追踪 方法 ， 研 究 者 可 以 建立 更 有 针对 性 的 、 针 


对 于 个 体 的 模型 ， 从 而 更 加 精准 地 识别 个 体 特征 3 


对 其 进行 更 有 效 的 治疗 (Wright & Woods, 


2020); 并 通过 探究 变量 在 被 试 间 水 平 的 差异 ， 检 验 单个 个 体 的 效应 是 否 能 够 推广 或 者 泛 化 


到 更 多 的 样本 中 (Hofmans et al., 2019). 


3 密集 追踪 数据 的 传统 分 析 方 法 


密集 追踪 数据 具有 典型 的 分 层 结构 ， 即 水 平 1( 时 间 点 水 平 ) 的 基本 单元 肉 套 于 水 平 2( 个 


体 水 平 ) 的 结构 ， 水平 2 的 结构 可 以 继续 藤 套 于 更 高 的 水 平 (如 ， 组 别 水 平 )， 因 此 传统 上 利用 


MLM 进行 分 析 ( 张 银 普 等 , 2016; Schwartz & Stone, 2007; Walls & Schafer, 2006)， 将 变量 的 变 


异 来 源 分 为 个 体内 和 个 体 间 等 多 个 水 平 (Walls & Schafer, 2006). 


3.1 方法 原理 


以 两 层 的 MLM 为 例 , 一般 来 说 ， 密 集 追 踪 数 


中 的 水 平 1 关注 个 体内 变量 在 不 同时 间 


点 之 间 的 关系 ， 而 水 平 2 关注 个 体 水 平 上 的 变量 间 


2015): 


Yit = Poi 十 B, Xi 


的 关系 。 对 于 水 平 1( 时 间 点 水 平 )， 以 存 


在 一 个 水 平 1 预测 变量 的 情况 为 例 , 个 体 i 在 测量 时 间 点 t 的 模型 公式 如 下 (Heck & Thomas, 


t + Eit (1) 


其 中 ,Xi 表示 个 体 ;在 第 上 个 测量 时 间 点 的 测量 ， 扩 ,和 ,分 别 表示 截 距 和 斜率 ,表示 误差 


项 (six~N(0,o?))。 男 外 ， 在 密集 追踪 数据 的 模型 中 ， 时 间 Tit 可 以 作为 预测 变量 加 入 到 水 平 1 


的 模型 中 。 


对 于 水 平 2( 个 体 水 平 )， 同 样 以 存在 一 个 水 平 2 预测 变量 的 情况 为 例 ， 其 模型 表达 式 如 


Boi = Yoo + HoWi + toi (2) 


Pii = Nort %1 Wit ti (3) 


其 中 ， ANB, ,分 别 为 水 平 1 中 的 随机 截 距 和 随机 斜率 ，Wi 为 水 平 2 IER, Ao ha 
分 别 为 水 平 2 的 预测 变量 Wi 对 BP,、pBi, 的 斜率 ， yo。、1 表 示 对 应 的 截 距 。 各 变量 在 水 平 1 的 
斜率 和 截 距 可 能 在 不 同 个 体 间 存在 变异 (随机 效应 )， 且 这 一 变异 可 能 可 以 被 水 平 2 的 预测 变 
量 解释 。 另 外 , 水平 1 和 水 平 2 的 模型 均 可 以 处 理 含 有 多 个 预测 变量 的 情况 。 根 据 实 际 研究 
中 密集 追踪 数据 的 具体 情况 ， 可 以 对 其 进行 更 多 层 的 建 模 ， 如 测量 时 间 点 嵌 套 于 个 体 ,个 体 
套 于 所 在 的 群体 等 。 

在 使 用 MLM 分 析 传 统 的 追踪 数据 时 ， 通 常 假设 不 同 个 体 在 水 平 1 的 残 差 方差 o? 均 相 
等 ， 但 在 密集 追踪 数据 中 ， 这 种 假设 很 可 能 会 带 来 一 些 问 题 (Jahng & Wood, 2017)， 如 参数 标 
准 误 的 有 偏 估计 等 (Hamaker & Grasman, 2015)。 密 集 追 踪 数 据 中 多 次 重复 测量 的 时 间 间 隔 较 
短 , 个 体 在 当前 时 间 点 的 测量 很 可 能 会 受到 先前 测量 的 影响 , 因此 通常 需要 假设 水 平 1 的 残 
差 带 有 自 回归 的 效应 (Hamaker & Grasman, 2015)。 以 考虑 一 阶 自 回 归 效 应 (Lag-1 


Autoregressive Effects) 为 例 ， 此 时 式 (1) 中 的 残 差 & 可 以 表示 为 (Hamaker & Grasman, 2015): 


ii 
t 
g 
it 


a 


El 


Eit = Déit_1 + eit (4) 
其 中 , 残 差 项 er 可 以 视 为 模型 中 没有 考虑 到 的 因素 的 影响 , 服从 正 态 分 布 , 即 ett~N(0，ae2)。 
DRIM i HÄ EIR, 即 该 个 体 在 当前 测量 时 间 点 的 测量 偏差 & 在 多 大 程度 上 受到 上 
一 测量 时 间 点 的 影响 。 自 回归 系数 GB 在 水 平 2 的 表达 式 为 (Hamaker & Grasman, 2015): 


®©; = Ø+ U2; (5) 
其 中 ，g 表 示 所 有 个 体 平均 的 自 回归 效应 ，uzi; 表 示 个 体 i 的 自 回 归 效 应 与 平均 自 回归 效应 
之 间 的 偏 移 ， 服 从 正 态 分 布 ，vzi~NM(0，c22)。 


© 3.2 传统 多 层 线 性 模型 的 不 足 


传统 的 MLM 可 以 解释 不 同 个 体 间 变量 随时 间 变 化 的 关系 曲线 存在 的 差异 (Walls & 
Schafer, 2006)， 但 是 这 种 模型 不 够 灵活 。 在 分 析 动 态 变化 效应 、 探 究 变 量 间 关 系 的 强度 和 方 
向 随时 间 变 化 的 趋势 时 ， 传 统 的 MLM 仍 存在 一 定 的 局 限 (Lietal.,2006)。 首 先 ， 当 研究 问题 
关注 的 是 个 体内 心理 过 程 的 动态 变化 ， 如 变异 性 和 不 稳定 性 时 , 传统 的 MLM 方法 在 分 析 个 
体内 发 展 模式 中 的 个 体 间 差异 时 可 能 会 产生 不 正确 的 参数 估计 (Jahng, 2008)， 且 不 能 解决 变 
量 随时 间 变 化 的 轨迹 在 个 体 间 存在 较 大 差异 的 问题 (Shiyko et al., 2012)。 其 次 ，MLM 通常 假 
设 随机 效应 服从 一 定 的 分 布 类 型 ， 如 正 态 分 布 ; 但 通常 情况 下 ， 这 种 假设 很 难得 到 满足 ， 因 


此 基于 这 些 假设 得 到 的 模型 参数 也 可 能 是 不 准确 的 (Piccirillo & Rodebaugh, 2019)。 针 对 这 一 


问题 , Ke 等 人 (2019) 指 出 , 将 贝 叶 斯 多 层 结构 方程 模型 (Bayesian Multilevel Structural Equation 


Modeling) 与 更 稳健 的 、 使 用 非 正 态 分 布 (如 Student's t 分 布 , Student's ¢ Distribution; 混合 正 态 
分 布 , Mixture of Normal Distributions; Tong & Zhang, 2012; Zhang, 2013; Zhang et al., 2013) 的 
贝 叶 斯 估计 方法 相 结合 ， 可 能 可 以 解决 随机 效应 不 符合 正 态 分 布 时 参数 估计 不 准确 的 问题 。 

另外 ， 在 MLM 中 ， 研 究 者 往往 需要 通过 中 心 化 来 分 解 组 间 水 平和 组 内 水 平 的 效应 
(Curran & Bauer, 2011)。 传 统 上 对 于 预测 变量 的 处 理 方法 包括 个 体 均值 中 心 化 、 总 均值 中 心 
化 等 ， 但 这 些 方法 的 适用 范围 较为 有 限 ， 在 一 些 情况 下 容易 产生 较 大 的 偏差 (Asparouhov & 
Muthén, 2019a; Curran & Bauer, 2011)。 如 ， 在 存在 情境 效应 2(Contextual Effects; Liidtke et al., 
2008) 或 自 回归 效应 (Nickel 1981)， 以 及 预测 变量 为 (两 个 及 以 上 类 别 的 ) 分 类 变量 且 存 在 非 
线性 连接 函数 (Link Function) 的 情况 下 ， 传 统 方法 均 会 产生 较 大 的 参数 估计 偏差 (Asparouhov 
& Muthén, 2019a)。 当 缺失 数据 是 随机 缺失 (Missing at Random, MAR) 而 不 是 完全 随机 缺失 
(Missing Completely at Random, MCAR) 时 ，MLM 中 的 样本 均值 无 法 体现 预测 变量 真实 的 均 
值 ， 传 统 的 中 心 化 处 理 同样 可 能 带 来 参数 估计 上 的 偏差 (Asparouhov & Muthén, 2019a)。 


4 密集 追 踊 数 据 分 析 方 法 的 新 进展 


为 了 克服 传统 方法 在 密集 追踪 数据 中 的 局 限 性 , 帮助 研究 者 更 好 地 挖掘 个 体内 部 动态 变 
化 的 过 程 和 作用 机 制 ， 下 面 将 主要 介绍 两 种 较 新 的 密集 追踪 数据 的 分 析 方 法 ， 包 括 
DSEM(Asparouhov et al., 2018) 和 结合 了 个 体 和 群体 信息 的 GIMME(Gates & Molenaar, 2012). 
根据 Hamaker 等 人 (2018) 的 界定 , 多 个 体 的 密集 追踪 数据 分 析 方法 可 以 分 为 自 上 而 下 的 方法 
(Top-Down Approach) 和 自 下 而 上 的 方法 (Bottom-Up Approach) 两 大 类 。 自 上 而 下 的 方法 通过 
较 高 水 平 的 模型 向 下 对 较 低 水 平 的 动态 效应 进行 解释 (如 MLM，DSEM 等 )， 这 里 的 高 低 次 
序 均 为 相对 的 而 不 是 绝对 的 高 低 ( 如 在 三 水 平 的 MLM 中 ,水 平 2 相对 于 水 平 1 是 较 高 水 平 ， 
但 相对 水 平 3 是 较 低 的 水 平 )。 而 自 下 而 上 的 方法 首先 对 每 个 个 体 的 数据 进行 分 析 ， 再 进 一 
步 向 上 搜寻 不 同 个 体 间 相 似 的 模式 (如 GIMME). 
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41 自 上 而 下 的 方法 : 动态 结构 方程 模型 


DSEM 承接 了 传统 MLM 自 上 而 下 的 思想 , 同时， 由 于 它 采 用 了 贝 叶 斯 的 框架 ， 因 此 相 


比 传统 MLM 更 为 灵活 (Asparouhov et al., 2018)。DSEM 结合 了 多 种 模型 方法 的 优势 ， 包 括 


MLM、 时 间 序 列 模型 (Time-Series Modeling)、 结 构 方程 模型 (Structural Equation Modeling, 


SEM) 以 及 时 变 效 应 模型 (Time-Varying Effects Modeling, TVEM)， 被 认为 是 用 于 分 析 密 集 追 


? 即 ， 在 控制 了 水 平 1 的 自 变量 Xi 的 情况 下 ， 水 平 2 的 平均 水 平 妨 仍 对 因 变 量 存在 额外 的 效应 。 


踪 数据 最 有 前 景 的 方法 之 一 (Zhou et al., 2020)。 它 强调 变量 在 时 间 尺 度 上 的 动态 变化 过 程 ， 
这 里 的 “动态 ”可 以 理解 为 变量 在 当前 时 刻 的 测量 会 受到 变量 自身 在 先前 时 间 点 的 影响 ( 自 
回归 效应 )， 以 及 受到 其 他 变量 在 先前 (交叉 滞后 效应 , Cross-lagged Effects) 以 及 当前 时 间 点 
(同时 效应 , Contemporaneous Effects) 的 测量 的 影响 。 另 外 ， 在 DSEM 中 ， 各 变量 的 变异 均 可 
通过 个 体 特异 和 时 间 特 异 的 变量 进行 解释 。 


4.1.1 方法 原理 


DSEM 可 以 认为 是 动态 因子 模型 (Dynamic Factor Model, DFM; Ferrer & Nesselroade, 2003; 
Molenaar, 1985; Molenaar et al., 1999; Zhang & Nesselroade, 2007) 在 多 个 个 体 情 况 下 的 扩展 
(Asparouhov et al., 2018)。 它 结合 多 层 模 型 的 思路 ,将 第 i 个 个 体 在 第 t 个 时 间 节 点 的 观测 变量 
Vie 分 解 为 个 体 i 的 特异 贡献 i;、 时 间 t 的 特异 贡献 73, 以 及 除去 个 体 和 时 间 的 特异 贡献 后 ， 
个 体 i 在 时 间 节 点 t 的 波动 变化 天 二 这 三 个 成 分 ， 即 Ye = Yit + You + Y3,¢(Asparouhov et al., 
2018)。 例 如 ， 对 被 试 群体 在 工作 日 中 感受 到 的 工作 压力 和 积极 情绪 水 平 进行 测量 ， 其 中 个 
体 i 在 时 间 节 点 t 的 工作 压力 水 平 7i,:， 将 同时 取决 于 个 体 i 在 工作 压力 上 的 一 般 基线 水 平 六;、 
时 间 节 点 t 对 个 体 工 作 压 力 的 一 般 效 应 (如 , 一 般 来 说 , 每 到 周一 时 个 体 的 工作 压力 都 会 增 
高 )， 以 及 除去 个 体 和 时 间 的 一 般 性 效应 后 ， 个 体 i 在 时 间 节 点 t 的 工作 压力 变化 ix( 如 公式 
(6)-(7) 所 示 ， 这 些 变 化 的 来 源 可 能 包括 : 该 个 体 在 该 时 间 节 点 的 积极 情绪 变化 ， 先 前 时 间 点 
的 积极 情绪 和 工作 压力 水 平 对 当前 时 间 点 的 影响 等 )。 

完整 的 DSEM 模型 的 框架 包括 1;、 妈 以 及 Yiit 这 三 个 成 分 ， 但 不 是 所 有 的 密集 追踪 

究 均 存在 时 间 特 异 的 效应 ， 因 此 不 是 所 有 的 密集 追踪 研究 均 适 合 采 用 最 完整 的 DSEM。 
Q 如 ， 在 密集 追踪 的 研究 中 ， 如 果 t = 1 仅 表示 个 体 的 第 一 个 观测 时 间 点 的 测量 ， 而 不 是 某 个 
特定 的 时 间 节 点 (比如 t = 1 都 表示 每 周 的 第 一 天 ) 时 ，t = 1 对 于 每 个 个 体 来 说 没有 一 致 的 含 
义 ， 因 此 数据 本 身 很 可 能 不 存在 时 间 特 异 的 效应 (Asparouhov et al., 2018). Asparouhov 等 人 
(2018) 指 出 , 研究 者 可 以 根据 实际 研究 的 需要 而 舍弃 ,这 一 成 分 , 只 考虑 个 体 水 平 的 随机 效 
应 ， 即 采用 两 层 的 DSEM(Two-Level DSEM); 甚至 进一步 舍弃 ;， 只 考虑 单 层 的 DSEM 模 
型 (Single-Level DSEM)。 两 层 和 单 层 的 DSEM 均 为 总 框架 下 的 特殊 的 情况 ， 因 此 ， 本 文 将 对 

总 的 DSEM 框架 下 的 方法 原理 进行 介绍 

首先 ， 与 DFM 类 似 ，DSEM 在 水 平 1 考虑 了 变量 在 时 间 节点 上 的 测量 可 能 受到 当前 及 


先前 工 次 测量 的 影响 ， 包 括 测量 模型 和 结构 方程 模型 (Asparouhov et al., 2018): 
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矩阵 。eLit、&iit 表 示 残 差 ， 均 值 均 为 0。 另 外 ，DSEM 框架 在 水 平 2 及 水 平 3 的 模型 公式 如 


下 (Asparouhov et al., 2018): 


Yoi = V2 + Aon, , + K2X21 + 2i (8) 
No = % + Bon, + T2X2i + zi (9) 
yat = V3 + Asnse + K3X34 + Ezt (10) 
R Nz, = 43 + Bang, +T3X3t + $3, (1 1) 


= EP, X PARR CEEZR A PE SE, (ELIS TARR, XG, 9 — 
组 时 间 特异 (在 不 同时 间 节 点 间 存在 差异 ， 但 不 随 个 体 变化 ) 的 变量 。n2; 和 ns 分 别 为 和 
凤 : 对 应 的 潜 变 量 。 总 的 来 说 ，DSEM 在 水 平 1 考虑 了 时 间 节 点 t 的 测量 与 当前 及 先前 次 测 
量 的 回归 关系 ， 在 水 平 2 考虑 个 体 特异 的 效应 ;而且 ， 与 传统 多 层 模型 不 同 的 是 ，DSEM 还 
= 在 水 平 3 考虑 了 时 间 特 异 的 效应 。 

= 具体 来 说 , DSEM 在 水 平 1 中 的 回归 系数 、 因 子 负载 和 截 距 项 等 参数 的 随机 效应 可 以 进 
© 步 在 水 平 2 和 水 平 3 进行 解释 ， 同 时 考虑 不 同 个 体 和 时 间 节点 的 特异 贡献 (Asparouhov et 
al.,2018)。 对 于 水 平 1 的 系数 s( 如 , s 可 以 指 代 回 归 系数 矩阵 Riz 中 的 元 素 、 因 子 负载 矩阵 Ai u 
中 的 元 素 、 截 距 项 ea ie 等 ), 其 随机 效应 可 以 在 水 平 2 和 水 平 3 被 分 解 为 sz; 和 sat 两 个 成 分 ， 
并 分 别 被 个 体 特异 和 时 间 特 异 的 变量 解释 (Asparouhov et al., 2018)。 如 ， 在 最 完整 的 DSEM 
框架 下 , 个 体积 极 情绪 的 自 回归 系数 的 随机 效应 将 被 分 解 为 个 休 和 时 间 两 个 方面 的 来 源 。 这 
种 随机 效应 在 水 平 2 可 能 受到 人 格 特质 (个 体 特异 的 变量 ) 的 影响 ， 同 时 在 水 平 3 受到 特定 工 
作 日 (时 间 特 异 的 变量 ) 的 影响 。 类 似 地 ， 对 于 模型 中 的 残 差 ， 其 方差 的 随机 效应 也 可 以 抽象 
地 表示 为 指数 函数 的 形式 = Exp(s2; + sst)， 从 而 可 以 分 别 利 用 个 体 和 时 间 特异 的 变量 对 


其 随机 效应 进行 解释 (Asparouhov et al., 2018). 


= 


McNeish 和 Hamaker(2019) 指 出 ， 相 比 于 MLM, DSEM 的 亮点 在 于 它 可 以 很 方便 地 考 
虑 变量 的 自 回归 效应 和 交叉 滞后 的 效应 。 标准 的 MLM 软件 也 可 以 创造 滞后 时 间 点 的 新 变量 
加 入 模型 ， 从 而 进行 和 DSEM 类 似 的 分 析 ; 但 当时 间 间 隔 不 等 或 存在 缺失 数据 时 , MLM 中 


的 缺失 数据 点 也 随 之 增多 ， 从 而 对 参数 估计 带 来 更 严重 的 问题 (McNeish & Hamaker, 2019). 
而 DSEM 是 利用 卡尔 曼 滤 波 (Kalman Filter; Harvey, 1989)， 基 于 上 一 时 间 点 的 观测 对 当前 缺 
失 时 间 点 的 数据 迭代 地 进行 预测 和 更 新 ， 不 需要 依赖 按 列 删 除 (Listwise Deletion)， 因 此 在 处 
理 缺 失 数据 方面 有 较 大 的 优势 。 另 外 ， 如 3.2 节 所 述 ， 在 MLM 中 ， 研 究 者 往往 需要 通过 中 
心 化 来 分 解 组 间 水 平和 组 内 水 平 的 效应 (Curran & Bauer, 2011)， 但 传统 的 显 变量 中 心 化 在 密 
追踪 数据 中 表现 不 佳 (Asparouhov & Muthén, 2019a)。DSEM 允许 对 潜 变 量 进行 建 模 ， 采 用 
潜 变 量 中 心 化 的 方法 ( 即 在 结构 方程 模型 部 分 自 变 量 X 被 分 为 XLit、X2i;、Xs.t 三 个 成 分 分 别 
进行 建 模 ) 来 排除 来 源 于 个 体 和 时 间 的 一 般 性 影响 (Asparouhov et al., 2018). 在 大 多 数 情况 下 ， 
这 种 潜 变量 中 心 化 的 方法 在 参数 估计 方面 的 表现 均 比 传统 显 变量 中 心 化 的 方法 更 优 , 能 很 好 
地 解决 传统 的 显 变量 中 心 化 带 来 的 问题 (Asparouhov & Muthén, 2019a)。 

需要 注意 的 是 ，DSEM 要 求 时 间 节 点 t 为 整数 ， 但 在 实际 数据 中 这 一 条 件 很 可 能 无 法 满 
足 (Asparouhov et al., 2018)。 例 如 ， 各 测量 节点 之 间 的 时 间 间 隔 可 能 参差 不 齐 (West, 2018), 
使 得 时 间 节 点 t 不 为 整数 。 为 了 解决 这 一 问题 DSEM 需要 建立 另 一 个 度量 网 络 ， 对 时 间 节 
点 进行 转换 和 处 理 ， 从 而 得 到 一 系列 近似 的 整数 值 (Asparouhov et al., 2018)。 例 如 ， 当 四 个 
测量 节点 的 间隔 分 别 为 0.5 天 、0.9 天 和 2.1 天 ， 可 以 选择 一 个 较 小 的 时 间 间 隔 6 = 0.5 作为 
网 格 的 最 小 单位 ， 将 四 个 测量 时 间 节 点 近似 转化 为 如 = 0, t= 1, = 2,6= 40 BF ts PEE 
观测 值 , DSEM 将 该 节点 视 为 缺失 值 处 理 (Asparouhov et al., 2018). Asparouhov 等 人 (2018) 指 
出 , 过 小 的 时 间 间 隔 6 会 导致 较 多 的 缺失 值 ， 而 间隔 偏 大 则 可 能 会 导致 多 个 时 间 节点 转化 为 
同一 个 整数 ， 因 此 需要 慎重 考虑 6 的 设置 。 对 于 如 何 选择 较为 合适 的 时 间 间 隔 5 值 ， 
Asparouhov 等 人 (2018) 提 供 了 较为 详细 的 建议 。 

一 般 来 说 ,， 贝 叶 斯 估计 不 需要 依赖 渐进 理论 , 在 小 样本 的 情况 下 相 比 频率 学 派 的 方法 具 
有 潜在 的 优势 (van de Schoot et al., 2014)， 但 这 种 优势 不 是 仅 通过 单纯 改变 估计 方法 便 能 得 
以 实现 的 (Depaoli & Clifton, 2015; van de Schoot et al., 2015)， 往 往 还 需要 研究 者 结合 实际 情 
况 对 参数 设置 合适 的 先 验 分 布 。 与 其 他 许多 贝 叶 斯 框架 下 的 模型 类 似 (Depaoli & van de 
Schoot, 2017; McNeish, 2017; Zondervan-Zwijnenburg etal., 2017), 在 样本 量 较 小 时 使 用 Mplus 
默认 的 无 信息 先 验 会 使 得 DSEM 更 倾向 于 产生 有 偏 的 估计 结果 (McNeish, 2019)。 因 此 ， 在 
样本 量 较 小 时 , 建议 研究 者 尽 可 能 结合 前 人 研究 和 其 他 先 验 知识 对 参数 设置 有 信息 的 先 验 分 
布 ; 在 缺乏 先 验 信息 的 情况 下 , 即使 是 根据 参数 取 值 范围 的 直观 判断 来 对 参数 设置 容许 范围 


限制 的 先 验 分 布 (Admissible-Range-Restricted Priors)， 也 可 以 较 大 程度 上 改善 DSEM 的 参数 


wun 


估计 结果 (McNeish, 2019). 


4.1.2 模型 拟 合 及 比较 


DSEM 可 以 利用 偏差 信息 准则 (Deviance Information Criterion, DIC) 作 为 模型 拟 合 与 比较 


的 指标 (Asparouhov et al., 2018). DIC 常用 于 贝 叶 斯 估计 中 的 模型 比较 (Spiegelhalter et al., 


2002), DIC 越 小 的 模型 


! 拟 合 越 好 。 然 而 ， 当 DSEM 中 参数 数量 较 多 、 模 型 过 于 复杂 时 ， 


要 较 多 迭代 次 数 才能 得 到 较为 稳定 的 DIC(Asparouhov et al., 2018)。 因 此 ，Asparouhov 等 人 


(2018) 建 议 , 研究 者 需要 通过 对 模型 设置 不 同 的 随机 种 子 并 多 次 运行 模型 ， 以 确认 DIC 是 否 


已 经 达到 稳定 的 取 值 。 


4.1.3 优势 及 不 足 


DSEM 在 密集 追踪 数据 的 分 析 上 存在 诸多 亮点 和 优势 。 首先 , DSEM 结合 多 种 模型 方法 
的 优势 ， 能 更 灵活 、 方 便 地 对 模型 中 多 个 变量 之 间 建 立 相 关 关 系 (Hamaker et al., 2018; Zhou 


etal., 2020)。 它 不 仅 局 限于 对 两 变量 的 建 模 分 析 ， 也 可 以 对 三 个 及 以 上 变量 进行 建 模 。 不 管 
是 显 变量 还 是 潜 变 量 、 是 否 是 在 同一 时 间 节 点 测量 , DSEM 中 各 个 变量 之 间 均 可 建立 回归 


径 (Asparouhov et al., 2018)。 其 次 ，DSEM 能 方便 地 将 密集 追踪 数据 中 呈现 的 特征 纳入 模型 


ms 


(如 ， 考 虑 数据 中 存在 的 动态 趋势 、 周 期 变化 等 )， 不 仅 可 以 考虑 变量 本 身 随时 间 的 变化 ， 还 
可 以 考虑 变量 间 关 系 随 时 间 的 动态 变化 (Asparouhov et al., 2018; Zhou et al., 2020)。 另 外 ， 


DSEM 对 模型 参数 的 估 
对 于 复杂 模型 的 处 理 更 


计 采 用 贝 叶 斯 模型 的 框架 , 相 比 于 传统 的 频率 学 派 方 法 ， 贝 叶 斯 方法 


有 优势 ( 张 沥 今 等 , 2019; Zhou et al., 2020). Song 和 Zhang(2014) 曾 将 


多 层 结构 方程 模型 与 DFM 相 结 合 ， 提 出 多 层 的 动态 因子 模型 (Multilevel Dynamic Factor 
Model, MDFM), 但 当 MDFM 中 存在 随机 效应 时 ,往往 会 导致 模型 无 法 估计 的 问题 ,而 DSEM 
对 于 测量 节点 和 随机 效应 的 设置 没有 数量 的 限制 ， 即 使 存在 较 多 的 测量 时 间 节 点 和 随机 效 


应 , 也 不 会 对 其 带 来 较 多 计算 量 上 的 负担 (Asparouhovetal.,2018)。 前 人 模拟 研究 表明 , DSEM 


在 随机 效应 的 均值 估计 


方面 表现 十 分 优秀 (Schultzberg & Muthén, 2018). 


然而 ，DSEM 也 存在 一 些 不 足 与 局 限 。 在 对 于 间隔 不 等 的 时 间 节 点 的 处 理 中 , 不同 时间 


间隔 5 的 选择 会 对 DSEM 水 平 1 参数 的 估计 结果 产生 一 定 的 影响 (Asparouhov et al., 2018; 


Asparouhov & Muthen, 2019b). Asparouhov 等 人 (2018) 在 DSEM 的 基础 上 做 出 微调 ， 提 出 了 


残 差 动态 结构 方程 模型 
型 中 体现 自 回 归 和 交叉 消 
影响 (Asparouhov et al. 


节点 的 间隔 不 一 的 情况 


(Residual Dynamic Structural Equation Model , RDSEM)， 仅 在 测量 模 

带 后 的 效应 , 因此 其 结构 方程 模型 部 分 不 会 受到 不 同时 间 间 隔 选择 的 
, 2018)。Asparouhov 和 Muthén (2019b) 的 模拟 研究 表明 ， 在 观测 时 间 
F, RDSEM 的 模型 结果 相对 DSEM 来 说 更 加 稳健 ， 几 乎 不 会 受到 时 


间 间 隔 的 影响 。 
总 的 来 说 ，DSEM 在 密集 追踪 数据 的 分 析 上 存在 诸多 亮点 和 优势 (Hamaker et al., 2018), 
但 也 仍 面临 许多 新 的 挑战 ， 亟 待 更 多 心理 学 、 统 计 学 等 领域 的 研究 者 共同 解决 。 


4.2 自 下 而 上 的 方法 : 组 迭代 多 模型 估计 


本 质 上 来 说 ，DSEM 和 传统 MLM 一 样 ， 将 被 试 视 为 整体 的 、 同 质 化 的 群体 ， 认 为 在 个 
体 间 存 在 相同 的 模型 路 径 ， 不 同 被 试 之 间 可 能 只 是 效应 大 小 有 所 不 同 (这 些 随机 效应 可 以 利 
用 个 体 特异 和 时 间 特 异 的 变量 进行 解释 ), 因此 是 一 种 基于 整体 的 方法 (Piccirillo & Rodebaugh, 
2019)。 然 而 ， 基 于 整体 和 基于 个 人 的 分 析 方 法 可 能 得 到 截然 不 同 的 结果 (Fisher et al., 2018; 
Turner & Hayes, 1996)， 特 别 是 实际 中 个 体 特异 性 较 强 的 情况 下 ， 针 对 整体 数据 构建 的 群体 
模型 可 能 不 能 解释 群体 内 部 的 任 一 个 体 (Miller & Van Horn, 2007; Molenaar & Campbell, 2009). 
即使 患 有 同一 疾病 ， 不 同 个 体 之 间 的 具体 情况 、 变 化 模式 等 也 存在 较 大 的 个 体 差 异 (Hayes et 
al., 2019)。 因 此 ， 近 年 来 越 来 越 多 学 者 倡导 研究 者 开展 基于 个 体 的 研究 和 统计 分 析 (Foster & 
Beltz, 2018; Hayes et al., 2019; Piccirillo & Rodebaugh, 2019; Wright & Woods, 2020)。 利 用 基于 
个 人 的 分 析 方 法 对 不 同时 间 、 不 同情 境 下 的 动态 变化 过 程 进行 建 模 和 分 析 , 可 以 得 到 更 精确 
的 、 更 具 针 对 性 的 结果 , 从 而 增加 治疗 和 干预 对 于 个 人 的 匹配 度 和 治疗 效用 (Hayes et al., 2019; 
Wright & Woods, 2020). 

纯粹 基于 个 人 的 分 析 方 法 充分 考虑 了 个 体 间 的 异 质 性 ， 但 同时 也 受到 了 一 些 质疑 。 如 ， 
仅仅 基于 个 人 建立 的 模型 可 能 不 够 可 靠 , 容易 受到 数据 中 噪音 的 影响 ， 从 而 得 到 过 拟 合 的 模 
型 (Gates & Molenaar, 2012)， 缺 乏 一 定 的 普 适 性 (Spencer & Schöner, 2003)。 由 此 ，Gates 和 


Molenaar (2012) 在 针对 个 人 的 模型 的 基础 上 ， 结 合群 体 共享 信息 的 优势 ， 提 出 GIMME 这 一 


算法 。GIMME 最 初 是 作为 一 种 针对 功能 性 磁 共 振 成 像 (Functional Magnetic Resonance 
Imaging, fMRI) 的 数据 分 析 方 法 提出 ， 但 近年 来 也 被 较 多 地 应 用 于 行为 和 心理 构 念 方面 密生 


追踪 数据 的 分 析 (Foster & Beltz, 2018; Gates et al., 2020)， 它 可 以 被 视 作 纯 粹 基于 个 体 的 方法 


由 


a 


和 基于 整体 的 方法 之 间 的 一 种 连接 (Beltz et al., 2016; Wright & Woods, 2020). 


4.2.1 方法 原理 


GIMME 是 一 种 动态 网 络 的 方法 ， 每 个 变量 均 被 视 为 网 络 中 的 一 个 节点 。 对 于 一 组 时 间 


序列 观测 变量 ni;(t)(wx1)，GIMME 的 公式 如 下 所 示 (Beltz & Gates, 2017): 


ni(t) = (A; + AZ) nit) + (Pii + DPn- 1) + EO (12) 


其 中 ，A(pxp) 表 示 相 同时 间 点 各 个 变量 节点 之 间 的 边 (Contemporaneous Edges), FJ LAE 


解 为 变量 间 的 同时 效应 ,89(pxp) 为 各 变量 节点 间 一 阶 的 滞后 效应 (First Order Lagged Edges), 
表示 在 时 间 点 上 时 ， 变 量 节 点 与 上 一 时 间 点 的 测量 之 间 的 关系 ， 可 能 同时 包括 自 回归 和 交叉 


灌 后 的 效应 。&,(t)(px1) 表 示 均 值 为 0 的 误差 项 。 


GIMME 的 个 体 模型 是 基于 统一 结构 方程 


模型 (Unified Structural Equation Modeling, USEM; Kim et al., 2007) 建 立 。 和 DFM 不 同 , USEM 


是 针对 显 变 量 ( 而 不 是 潜 变 量 ) 进 行 建 模 (Kim etal., 2007)。 在 USEM 的 基础 上 ，GIMME 还 区 


分 了 个 体 和 群体 的 效应 ， 下 标 i 表示 该 系数 为 第 i 


个 被 试 独 有 ， 而 上 标 g 表示 群体 模型 中 的 


效应 。 因 此 ， 在 各 个 个 体 的 模型 中 ， 将 同时 包括 个 体 特异 以 及 群体 共享 的 路 径 。 简 单 来 说 ， 


GIMME 的 建 模 分 析 步 又 大 致 包括 : (1) 数 据 预 处 型 


Rs (2) 基 于 大 多 数 个 体 模型 中 存在 的 路 径 ， 


自 下 而 上 地 建立 群体 的 模型 ， (3) 以 群体 模型 为 初始 模型 建立 个 体 的 模型 (Beltz & Gates, 2017; 


Gates & Molenaar, 2012)。 


GIMME 假设 认为 ， 残 差 项 E(t) 与 时 间 相 互 独立 ， 不 考虑 变量 随时 间 发 展 的 线性 趋势 


(Beltz & Gates, 2017)。 因 此 ， 研 究 者 需要 将 原始 数 


据 对 时 间 这 一 变量 进行 回归 分 析 ， 去 除 时 


间 对 变量 的 线性 影响 后 ， 再 将 得 到 的 残 差 项 纳入 GIMME 中 进行 正式 的 分 析 (Beltz & Gates, 


2017)。 另 外 ， 研 究 者 也 可 以 将 时 间作 为 外 源 变 量 (Exogenous Variable) 纳 入 模型 ， 从 而 去 除 各 


变量 随时 间 变 化 的 趋势 (Woods et al., 2020)。 


在 群体 模型 的 建立 中 ,GIMME 将 初始 模型 作为 起 点 ,根据 修正 指数 (Modification Indices) 


a 


对 个 体 模型 逐步 迭代 地 添加 回归 路 径 , 只 有 显著 优 


化 大 多 数 个 体 模型 拟 合 的 路 径 才能 被 纳入 


群体 模型 中 (Gates & Molenaar, 2012)。 接 着 ，GIMME 以 群体 模型 作为 迭代 的 初始 模型 ， 对 
各 个 个 体 模型 逐步 添加 显著 优化 其 模型 拟 合 的 路 径 ， 得 到 各 个 个 体 的 模型 。 关 于 GIMME 建 


立 群 体 和 个 体 模 型 的 具体 算法 可 见 Gates 和 Molenaar(2012) 的 文章 。 


一 般 来 说 ，GIMME 的 一 阶 延迟 模型 (包括 同时 效应 和 灌 后 一 个 时 间 点 的 效应 ) 已 经 能 较 


好 地 解释 密集 追踪 数据 中 的 时 间 依 赖 性 , 但 还 是 有 


(APosteriori Model Validatiom)， 以 查看 是 否 需要 考 


必要 对 每 个 个 体 的 模型 进行 后 验 模型 验证 


虑 更 高 阶 的 模型 (Beltz & Gates, 2017; Beltz 


& Molenaar, 2015)。 这 种 后 验 模型 的 验证 可 以 通过 检查 残 差 项 &.(t) 是 否 为 白 噪音 来 实现 (Beltz 


& Gates, 2017)。Beltz 和 Molenaar(2015) 详 细 介 绍 了 


并 提供 了 相应 的 LISREL 代码 。 


GIMME 中 后 验 模型 验证 的 方法 和 原理 ， 


在 GIMME 的 基础 上 ， 后 续 较 多 研究 也 提出 了 其 拓展 模型 ， 如 Gates 等 人 (2017) 提 出 的 


考虑 子 群 组 的 GIMME(Subgrouping Group Iterative Multiple Model Estimation, S-GIMME), 


Beltz 和 Molenaar(2016) 提 出 的 多 解决 方案 的 GIMME(GIMME for Multiple Solutions, GIMME- 


MS), LAR Gates 等 人 (2020) 提 出 的 基于 潜 变 量 进行 建 模 的 GIMME(Latent Variable GIMME, 


LV-GIMME). 在 GIMME 的 基础 上 , LV-GIMME 利用 带 有 工具 变量 (Instrumental Variables) 的 


两 阶段 最 小 二 乘法 估计 (Two-Stage Least Squares Estimation, MITV-2SLS)， 同 时 允许 测量 模型 
和 结构 模型 部 分 在 不 同 个 体 间 存在 差异 。 模拟 研究 的 结果 表明 , 相 比 于 对 条 目 总 分 直接 进行 


建 模 ，LV-GIMME 能 更 有 效 地 检测 出 模型 中 存在 的 效应 (Gates et al., 2020). 


4.2.2 优势 及 不 足 


总 的 来 说 ，GIMME 在 密集 追踪 数据 的 分 析 上 具有 其 独特 的 优势 。 与 基于 整体 的 分 析 方 
法 相 比 ，GIMME 的 优势 主要 体现 在 : (1) 对 于 群体 模型 的 建立 是 一 个 自 下 而 上 的 过 程 ， 其 路 
径 的 纳入 与 否 取 决 于 该 路 径 能 否 显著 改善 一 定 比例 的 单个 个 体 模型 (Gates & Molenaar, 2012), 
可 以 尽 可 能 减少 个 体 特异 信息 的 损失 , 更 适用 于 探究 变量 间 共 变 模式 的 个 体 差异 (Beltz et al., 
2016); (2) 人 允许 各 个 个 体 的 模型 结构 之 间 存 在 异 质 性 ， 并 且 能 提供 任 一 个 体 相 对 于 整个 群体 
的 差异 方面 的 信息 ,从 而 更 好 地 理解 个 体 模型 的 特征 和 模式 , 这 对 于 临床 上 针对 个 人 的 精准 
S 治疗 来 说 具有 重要 的 意义 (Wright & Woods, 2020)，(3) 对 于 被 试 数量 的 要 求 较 少 ， 在 只 有 10 
名 被 试 的 情况 下 就 已 经 可 以 较 好 地 对 群体 模型 的 路 径 进 行 估计 (Gates & Molenaar, 2012; 
Henry et al., 2019)。 而 与 纯粹 基于 个 体 的 方法 相 比 ，GIMME 也 有 具有 其 独特 的 优势 : 它 在 考 
ql 虑 个 体 特 异性 的 同时 大 大 提高 了 结果 的 普 适 性 ， 一 定 程度 上 降低 了 模型 过 拟 合 的 风险 (Beltz 
多 Gates,2017)。 另 外 ，GIMME 在 软件 操作 上 较为 简便 ， 且 能 自动 输出 每 个 个 体 模 型 以 及 群 
© 体 模 型 的 各 项 拟 合 指标 、 路 径 系数 等 ， 为 研究 者 基于 个 人 的 精细 化 研究 提供 了 丰富 的 信息 。 
= GIMME 还 可 以 较为 直观 地 对 结果 进行 呈现 , 更 易于 理解 与 解释 ， 适 用 于 动态 系统 中 存在 较 
< 多 相互 影响 的 变量 的 情况 (Beltz et al., 2016). GIMME 不 仅 可 以 帮助 研究 者 从 宏观 的 视角 理 
解 变量 间 的 关系 , 还 可 以 输出 每 个 个 体 的 模型 参数 估计 结果 ,大 大 增加 治疗 效用 和 干预 指导 
的 精确 性 。 
虽然 GIMME 在 对 密集 追踪 数据 的 分 析 上 存在 诸多 优势 ， 但 它 同 样 存在 一 些 局 限 。 首 
先 ，GIMME 本 质 上 是 一 种 数据 驱动 的 建 模 方法 ， 可 能 出 现 模 型 过 拟 合 (Overfitting) 等 问题 
(Beltz & Gates, 2017)。 其 次 ，GIMME 主要 考虑 变量 间 共 变 (Covary) 的 关系 及 机 制 ， 无 法 将 几 
乎 不 随时 间 变 化 的 变量 纳入 模型 进行 分 析 (Beltz & Gates, 2017)。 为 了 探究 个 体 水 平 上 不 随时 
间 变 化 的 协 变量 的 作用 ,研究 者 需要 额外 对 个 体 模型 中 得 到 的 参数 估计 进行 进一步 分 析 。 如 ， 
Kim 等 人 (2007) 利 用 一 般 线 性 模型 (General Linear Model), 探究 个 体 层 面 的 变量 (性 别 、 年 龄 
智商 及 教育 ) 对 USEM 中 得 到 的 个 体 模型 路 径 的 影响 。 又 如 ，Beltz 等 人 (2013) 利 用 t 检验， 


以 


比较 USEM 得 到 的 个 体 模型 路 径 中 ， 各 变量 间 达 到 显著 的 回归 路 径 的 数量 是 否 存在 显著 的 
性 别 差异 。 另 外 ，GIMME 假设 认为 变量 间 的 效应 相对 时 间 静 止 ， 不 考虑 数据 中 可 能 存在 的 
时 变性 质 (Time-Varying Nature; Beltz & Gates, 2017; Foster & Beltz, 2018)， 但 实际 研究 中 可 能 
出 现 变 量 间 关系 随时 间 动 态 变 化 的 情况 。 最 后 ，GIMME 要 求 各 个 测量 时 间 节 点 之 间 必 须 间 
阳 相 等 ,但 在 许多 密集 追踪 的 研究 中 ,观测 节点 的 时 间 间 隔 常 常 是 不 等 的 (Foster & Beltz, 2018; 
Kleiman et al., 2019)。 另 外 ， 虽 然 GIMME 对 于 个 体 数量 上 几乎 没有 限制 ， 但 它 要 求 对 于 每 


个 个 体 进 行 较 多 时 间 节 点 的 测量 。 


E 


4.3 实证 演示 及 方法 比较 
DSEM 已 经 可 以 在 Mplus 上 实现 ， 且 代码 语句 较为 简单 ， 这 为 研究 者 们 利用 DSEM 进 
行 密集 追踪 数据 的 动态 分 析 提 供 了 很 大 的 便利 (McNeish & Hamaker, 2019)， 近 年 来 涌现 出 较 


多 关于 DSEM 的 应 用 研究 (Armstrong et al., 2019; Hamaker et al., 2018; Racine et al., 2019; 


Zhang & Zheng, 2019). T GIMME 已 有 与 其 同名 的 及 (gimme; Lane etal,2020)， 也 可 以 利 


] LISREL 以 及 MATLAB 中 的 工具 箱 进行 建 模 ， 近 年 来 在 应 用 领域 也 出 现 较 多 研究 利用 


GIMME (Beltz & Gates, 2017; Bouwmans et al., 2018) 或 GIMME 的 扩展 模型 进行 密集 追踪 数 


据 的 分 析 ， 如 考虑 子 群 体 模 式 的 S-GIMME (Lane et al., 2019)， 以 及 可 以 同时 输出 多 个 群体 


模型 方案 的 GIMME-MS(Beltz et al., 2016). 


为 了 进一步 加 深 读 者 对 于 DSEM 和 GIMME 实际 建 模 过 程 的 理解 、 突 出 二 者 相 比 传统 
MLM 的 建 模 优势 ， 本 节 将 利用 Mplus 8.4 和 及 包 gimme, Xf McNeish 和 Hamaker(2019) 提 
供 的 公开 数据 分 别 进行 MLM, DSEM 和 GIMME 的 建 模 分 析 ( 具 体 代码 见 附录 )。 简单 来 说 ， 
McNeish 和 Hamaker(2019) 对 100 名 被 试 在 50 个 时 间 点 的 抑郁 水 平 (Depression) 和 吸烟 欲望 
(Urge to Smoke) 进 行 密集 追踪 ， 利 用 DSEM 一 系列 不 同形 式 的 模型 探讨 抑郁 水 平和 吸烟 欲 
望 之 间 的 动态 关系 ， 以 及 个 体 特异 的 变量 (工作 压力 和 家 庭 压力 ) 在 其 中 所 起 的 作用 。 


4.3.1 传统 多 层 线 性 模型 的 实证 演示 


McNeish 和 Hamaker(2019) 的 数据 中 同时 包括 个 体 层面 和 时 间 点 层面 的 变量 测量 ， 且 因 
变量 吸烟 欲望 的 组 内 相关 系数 (Intraclass Correlation, ICC) 较 大 (ICC = 0.387)， 因 此 可 以 利用 
MLM 进行 数据 分 析 。 假 设 吸烟 欲望 的 基线 水 平 、 随 时 间 变 化 的 趋势 以 及 抑郁 对 吸烟 欲望 的 
影响 在 不 同 个 体 间 存在 较 大 的 差异 ， 因 此 本 例 将 建立 带 随机 截 距 和 随机 斜率 的 MLM 模型 
(对 应 Mplus 代码 见 附录 A)。 为 了 更 好 地 解释 结果 ， 我 们 对 水 平 1 的 预测 变量 (抑郁 水 平 ) 进 
行 个 体 均 值 中 心 化 。 根 据 3.1 节 中 对 于 MLM 方法 原理 的 介绍 ， 在 水 平 1， 个 体 在 测量 时 


= 
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ELR t HO MAAR ERY, AY URRA: 


表示 测 


Yit = By; + B,(Xie — Xi) + PoiTie + Sie 


量 时 间 。 我 们 假设 水 平 1 各 项 斜率 和 截 距 在 不 同 个 体 间 存在 变异 , 且 这 一 变异 可 


以 被 个 体 本 身 来 自 工作 和 家 庭 的 压力 解释 。 因 此 ， 在 水 平 2，MLM 的 模型 表达 式 如 下 


其 


Boi = Yoo t HoW rii + HoW zgi + Uoi 


Pii = or F AW ries + Yor ži t i 


Boi = Moo + HaW rii + ha Wg 十 U2 


(13) 


其 中 ，Xiz 和 总 分 别 表 示 个 体 i 在 第 上 个 测量 时 间 点 的 抑郁 水 平 以 及 个 体 ; 的 平均 抑郁 水 平 ，T 


能 可 


(14) 
(15) 
(16) 


P, Woy, Wg SAKE 2 的 预测 变量 ， 分 别 表示 个 体 的 工作 压力 和 家 庭 压力 。 


MLM 对 应 的 模型 图 如 图 1 所 示 : 


ve: 圆圈 


> 


家 庭 压 力 ， 


1 抑郁 水 平 对 吸烟 欲望 的 影响 


内 变量 为 潜 变 量 ， 方 框 内 变量 为 显 变量 。 下 标 i 表示 第 i 个 个 体 的 测量 ，t 表示 第 上 个 时 间 点 的 观测 数据 。 


Mplus 8.4 默认 对 MLM 采用 稳健 极 大 似 然 估计 (Robust Maximum Likelihood Estimator, 


MLR), 
化 的 结 
更 强烈 


因此 仅 允 许 输 出 非 标准 化 的 结果 (Muthén & Muthén, 1998-2017)。 如 果 需 要 输出 
R, 可 以 采用 贝 叶 斯 估计 的 方法 。 结果 显示 ， 工作 压 力 和 家 庭 压力 越 高 的 个 体 往 
的 吸烟 欲望 (y 6=0.458, 7,=0.286, ps<0.001) 个 体 的 吸烟 欲望 随时 间 的 推移 而 逐 


加 (%。=0.157, p<0.001)， 但 这 种 随时 间 变 化 的 趋势 在 个 体 间 的 随机 效应 不 能 被 个 体 本 身 
作 压 力 或 家 庭 压力 解释 (ps>0.053)。 抑 郁 水 平 越 高 时 ， 个 体 的 吸烟 欲望 也 越 强烈 (7o; =0.741， 
p<0.001)， 且 个 体 的 工作 压力 和 家 庭 压力 越 大 ， 该 个 体 抑郁 水 平 对 吸烟 欲望 的 影响 也 


(X41=0. 270, p=0.001; 7,,=0.334, p<0.001)。 总 之 ， 个 体 本 身 的 工作 压力 和 家 庭 压力 会 提高 个 


体 吸 烟 


欲望 的 基线 水 平 ， 且 增加 其 抑郁 水 平 对 吸烟 欲望 的 影响 。 


标准 
往 有 
浙 增 


的 工 


越 大 


4.3.2 动态 结构 方程 模型 的 实证 演示 


与 MLM 相 比 ,， DSEM 可 以 方便 地 对 密集 追踪 数据 中 的 动态 效应 进行 建 模 分 析 。 本 例 仅 
考虑 了 个 体 特异 的 预测 变量 , 因此 采用 两 层 的 DSEM.. McNeish 和 Hamaker(2019) 的 数据 中 ， 
吸烟 欲望 整体 呈现 出 随时 间 线 性 增长 的 趋势 ， 因 此 需要 对 数据 进行 去 时 间 趋 势 (Detrending)。 
Asparouhov (2018) 的 模拟 研究 表明 , DSEM 在 同时 进行 去 时 间 趋 势 和 自 回归 效应 的 分 析 时 会 
对 方差 的 估计 产生 偏差 ， 因 此 建议 采用 DSEM 的 变 式 一 一 RDSEM 进行 分 析 。RDSEM 与 
DSEM 的 区 别 主要 在 于 RDSEM 分 别 对 去 时 间 趋 势 与 自 回 归 和 交叉 滞后 效应 进行 建 模 ， 即 
仅 在 测量 模型 部 分 而 不 是 结构 方程 模型 的 部 分 考虑 自 回 归 和 交叉 滞后 效应 。RDSEM 是 对 
DSEM 微调 后 的 提出 的 一 种 变 式 (Asparouhov et al., 2018)， 因 此 对 RDSEM 的 介绍 同样 可 以 


较 好 地 展示 DSEM 的 分 析 过 程 及 思路 。 


根据 理论 假设 ， 在 水 平 1， 吸 烟 和 欲望 和 抑郁 水 平 存在 自 回归 (@w ye gas Phonan) AE 


又 滞后 的 效应 , 即 前 一 个 时 间 点 的 抑郁 水 平 会 对 后 一 时 间 点 的 吸烟 欲望 产生 影响 (gw in a) 


而 前 一 时 间 点 的 吸烟 欲望 也 可 能 会 对 后 一 时 间 点 的 抑郁 水 平 产生 作用 (Byyiy_w)。 对 应 模型 


公式 为 : 
(w) 
al 7 Hesse 吸烟 站 (17) 
Yor ie) sana (w) 
抑郁 ,让 Y ttit 
(w) 
Ymie| Fy 《吸烟 让 Te 
YO) | [Bg-a oe 6 抑郁 it (18) 
抑郁 , 认 : g 
F ar B P Di ga k n y% | (19) 
Cigi) [Panes Pin- init- [Oana ie 


与 4.1.1 节 中 DSEM 的 公式 相 比 ，RDSEM 仅 在 测量 模型 部 分 考虑 残 差 项 5 的 自 回 归 和 


交叉 滞后 的 效应 。 令 Neg = Var (Spa) Tgi = Var 5 知名 站)， 则 在 水 平 2， 个 体 的 吸烟 


和 欲望 和 抑郁 的 平均 水 平 @ 豚 烟 Lyne EVAR Pye gsi Savina Bi- 
Da-mi Bam- Pygmy aya)? 2S BAe OT A yy yc? Magra, PI HESS BIT TE 


压力 和 家 庭 压力 的 影响 。 在 考虑 残 差 项 方差 的 随机 效应 时 ， 为 了 保证 残 差 项 方差 为 正 数 ， 


DSEM 将 取 对 数 之 后 的 残 差 项 方差 rg wg 和 cjgyfy i 纳入 水 平 2 的 模型 。 因 此 , 水平 2 模型 的 公 


式 为 : 


202009.00002v2 
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其 中 ， 公 式 (17) 对 应 表示 变量 变异 来 源 的 分 解 ， 公 式 (18)-(19) 对 应 水 平 1 的 模型 ， 公 式 
(20) 对 应 水 平 2 的 模型 (对 应 Mplus 代码 见 附录 A)。McNeish 和 Hamaker(2019) 的 数据 仅 提 供 
了 个 体 特异 预测 变量 的 数据 ， 因 此 本 例 仅 假设 水 平 1 的 各 项 系数 在 水 平 2 受到 个 体 特异 预 
测 变 量 的 影响 ; 但 如 果 数 据 信息 充足 , 我 们 也 可 以 类 似 地 在 水 平 3 考虑 时 间 特 异 的 预测 变量 


(如 ， 特 定 的 工作 日 ) 产 生 的 影响 ， 如 可 以 将 残 差 项 方差 ri 的 随机 效应 抽象 地 表示 为 


log (pass) = szit S3t， 从 而 将 其 分 解 为 个 体 和 时 间 特 异 的 随机 效应 。 在 本 例 中 , 对 于 第 


i 个 个 体 ， 其 模型 图 如 图 2 所 示 : 


变量 分 解 


吸烟 -吸烟 i 


多 抑郁 -吸烟 


人 Bl 吸烟 -抑郁 .i 


人 抑郁- 抑郁: 


2 吸烟 欲望 和 抑郁 水 平 的 动态 关系 


注 : 圆圈 内 变量 为 潜 变 量 ， 方 框 内 变量 为 显 变量 。(W) 表 示 个 体内 的 估计 ,“ 吸 烟 ” 和 “抑郁 ”分 别 表示 吸烟 欲望 和 抑郁 
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的 吸烟 欲望 对 当前 抑郁 水 平 、 上 一 时 间 点 的 抑郁 水 平 对 当前 吸烟 欲望 的 交叉 滞后 效应 。 UER Syn PS gga 


表示 吸烟 欲望 和 抑郁 水 平 的 残 差 项 。zaili 和 maniri 分 别 为 6uiz 和 afi 的 方差 。 为 了 更 清晰 地 呈现 模型 假设 ,模型 图 中 


未 呈现 与 去 时 间 趋 势 有 关 的 路 径 。 


DSEM 和 RDSEM 均 可 以 同时 输出 标准 化 和 非 标准 化 的 结果 ， 且 标准 化 结果 还 在 水 平 1 


输出 了 不 同 个 体 


Over Clusters). 


的 是 Mplus 输 


章 。 水 平 1 标准 化 系数 在 不 同 个 体 中 


(Dy sj —ygyg=0-177, CI=[0.142, 0.211])， 前 一 时 间 点 的 吸 类 


0.738])。 前 一 时 间 点 的 抑郁 水 平 对 后 一 时 间 点 吸烟 欲望 的 影 


为 了 更 


h 各 项 标准 化 系数 的 平均 结果 (Within-Level Standardized Estimates Averaged 
好 地 对 吸烟 欲望 和 抑郁 水 平 之 间 的 交叉 滞后 效应 进行 解读 ， 我 们 采用 
出 的 标准 化 结果 ， 其 标准 化 的 具体 过 程 及 原理 见 Schuurman 等 人 (2016) 的 文 
的 平均 结果 显示 ， 吸 烟 欲 望 的 自 回归 效应 达到 显著 


央 欲 望 可 以 显著 预测 后 一 时 间 点 的 吸 
烟 欲 望 。 男 外 ， 被 试 的 吸烟 欲望 整体 随时 间 的 推移 而 逐渐 增加 (By jy_wy=0. 721, CI=[0.704, 


啊 没 有 达到 显著 , 但 前 一 时 间 点 


的 吸烟 欲望 可 以 降低 后 一 时 间 点 的 抑郁 水 平 (By 各 _ 吸 灿 =-0.069, CI=[-0.102, -0.038])。 在 水 平 


2， 工 作 压 力 和 家 庭 压力 越 大 的 个 体 ， 其 吸烟 欲望 越 高 (yy wyw=0.362, CI=[0.190, 0.515]; 


Yzma =0-220, CI=[0.055, 0.371]】)， 吸 烟 欲 望 的 自 回归 效应 越 强 (yT 作 吸烟 -吸烟 =0.423， 


CI=[0.220, 0.594]; y 闪 庭 吸烟 -吸烟 ~0.212，CI=[0.016, 0.388]), MH 


(1 4,108( ty) 


=0.355, CI=[0.227, 0.474]; Y% 


We l0E( tg) OL CI=[0.076, 0.334])， 因 此 


欲望 的 残 差 项 方差 越 大 


吸烟 


欲望 的 残 差 的 动态 波动 幅度 也 越 大 。RDSEM 的 结果 进一步 表明 , 个体 的 吸烟 欲望 存在 较 强 


的 惯性 ( 自 回归 效应 )， 且 个 体 的 工作 压力 和 家 庭 压力 越 大 ，:; 


MLM 的 分 析 中 发 现 抑郁 水 平 对 吸烟 欲望 存在 正 向 预测 的 作用 ， 而 RDSEM 的 结果 体现 了 相 


反方 向 的 另 一 种 可 能 ， 


中 ， 变 量 间 的 相互 作用 很 可 能 存在 复杂 的 双向 反馈 ， 而 DSEM RÜ 
可 能 可 以 为 研究 者 带 来 更 多 新 的 研究 思路 和 启发 。 


的 作用 路 径 ， 这 


即 前 一 时 间 点 的 吸烟 欲望 可 以 降低 后 一 时 间 点 的 抑郁 水 平 。 现实 生活 


这 种 惯性 也 越 强 。 我 们 在 进行 


4.3.3 组 迭代 多 模型 估计 的 实证 演示 


本 节 将 利 月 


H gimme 这 一 R 包 进 行 GIMME 的 建 模 分 析 。 在 对 数据 进行 分 析 之 前 ， 


F 在 同一 模型 中 建立 双向 


首先 


需要 去 除数 据 中 随时 间 线 性 变化 的 趋势 ,研究 者 可 以 将 原始 数据 对 时 间 进 行 回归 后 得 到 的 残 


差 项 纳入 GIMME 中 进 


行 正式 的 分 析 , 也 可 以 直接 将 时 间作 为 外 源 变量 纳入 模型 ， 从 而 去 除 


各 变量 随时 间 变 化 的 趋势 (Woods et al.,2020)。 为 了 更 直观 地 呈现 去 时 间 趋 势 的 效果 ， 本 例 采 


前 一 种 去 时 间 趋 势 的 方法 。 以 第 一 个 被 试 的 吸烟 欲望 为 例 ， 图 3a 为 吸烟 欲望 的 原始 测量 
(URGE) 随 时 间 变 化 的 总 体 趋势 ， 随 着 时 间 的 推移 ， 吸 烟 欲 望 越 强烈 。 通 过 建立 时 间 对 吸烟 
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欲望 的 回归 模型 得 到 的 残 差 来 去 除数 据 中 随时 间 线 性 变化 的 趋势 后 ， 其 结果 如 图 3b 所 示 。 
可 以 看 出 ， 此 时 吸烟 欲望 在 不 同时 间 点 的 测量 呈 平 稳 的 趋势 。 


= 


URGE 


r 4 5 0 10 20 30 
TIME TIME 


3 吸烟 欲望 随时 间 发 生 的 变化 


注 : 图 3a 为 吸烟 欲望 的 原始 测量 随时 间 变 化 的 总 体 趋势 ， 图 3b 表示 去 除数 据 中 随时 间 线 性 变化 的 趋势 之 后 的 结果 。 


去 除 各 变量 受到 的 来 自 时 间 的 影响 后 ， 将 其 纳入 GIMME 进行 分 析 。 对 于 每 个 被 试 i, 
其 吸烟 欲望 和 抑郁 水 平 的 时 间 序 列 观测 变量 可 以 表示 为 : 


| _ | g Hed : g bo = st 天 | 
be (4: j Ai ) Naep,i (t) g (er i Pii Naep,i (t T 1) Sdep, i(t) (21) 


其 中 ，7uroei(b 和 7aepi(b 分 别 表示 被 试 在 时 间 点 ! 的 吸烟 欲望 及 抑 帮 水平。4(2x2) 表 
示 吸 烟 欲望 和 抑郁 水 平 的 同时 效应 矩阵，S@(C2x2) 为 吸烟 欲望 与 抑郁 水 平 之 间 的 一 阶 滞后 效 
应 (First Order Lagged Edges)， 包 括 自 回归 和 交叉 灌 后 效应 。 下 标 i 表示 该 系数 为 第 i 个 被 试 
独 有 ， 而 上 标 g 表示 为 群体 模型 中 的 效应 。 在 利用 GIMME 进行 数据 分 析 时 ,研究 者 可 以 利 
用 已 有 的 经 验 对 群体 模型 中 可 能 存在 的 路 径 做 出 假设 ， 并 将 其 作为 初始 迭代 的 模型 (如 ， 在 
gimmeSEM 函数 中 声明 paths=urge~deplag， 则 表示 初始 模型 中 存在 前 一 时 间 点 的 抑郁 水 平 
指向 当前 吸烟 欲望 的 回归 路 径 )。 若 理论 依据 不 足 ， 也 可 以 直接 从 空 模型 出 发 ， 以 数据 为 驱 
动 进行 探索 式 的 建 模 (Beltz & Gates, 2017; Gates & Molenaar, 2012)。 本 文 将 利用 数据 驱动 、 
探索 式 的 方式 进行 演示 。 

在 准备 好 数据 并 运行 函数 后 ，gimme 将 初始 模型 作为 起 点 ， 自 下 而 上 地 建立 群体 模型 。 
具体 来 说 ,gimme 首先 根据 修正 指数 对 每 个 个 体 模型 逐步 迭代 地 添加 回归 路 径 , 再 以 个 体 模 
型 中 的 路 径 为 基础 建立 群体 模型 。 为 了 找到 在 大 多 数 个 体 模 型 中 较为 稳定 存在 的 路 径 和 效 
应 ， 该 路 径 只 有 能 显著 优化 大 多 数 个 体 模型 的 拟 合 时 (大 多 数 个 体 模 型 中 ， 该 路 径 对 应 的 MI 
在 0.01 的 水 平 上 达到 显著 )， 才 能 被 纳入 群体 模型 中 (Gates & Molenaar, 2012)。 本 文 按照 默 
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认 的 设置 ， 将 75% 设 为 可 以 被 纳入 群体 模型 的 临界 值 ， 研 究 者 也 可 以 根据 具体 需要 进行 修 

改 。 基 于 群体 模型 ，gimme 再 对 个 体 模型 逐个 进行 迭代 路 径 搜索 。 在 完成 所 有 模型 的 路 径 搜 

RJ gimme 将 输出 各 个 模型 的 路 径 图 及 详细 的 统计 结果 。 图 4 即 gimme 输出 的 群体 模型 

图 。 群 体 模型 结果 显示 ， 超 过 75% 的 个 体 模型 中 ， 吸烟 欲望 和 抑郁 水 平均 存在 自 回 归 的 效 

应 ( B urge>urgelag=.196, Ê dep3deplag=-.049)。 而 其 他 路 径 均 没有 在 超过 75% 的 个 体 模型 中 出 现 ， 

这 一 结果 也 体现 了 变量 间 关 系 的 个 体 差异 。 
‘ 


图 4 群体 模型 


注 : urge: 吸烟 欲望 ，dep: 抑郁 水 平 。 黑 线 表 示 群 体 模 型 中 的 路 径 ， 灰 线 表 示 个 体 模型 中 存 如 


E 的 路 径 ， 实 线 代 表 同 时 效 


应 ， 虚 线 表 示 滞 后 效应 (包括 自 回归 和 交叉 滞后 效应 )。 线 条 粗细 表示 模型 中 存在 该 路 径 的 个 体 数 量 的 多 少 。 


图 5 为 随机 选取 的 个 体 模型 示例 。 总 体 来 看 , 不 同 个 体 在 吸 


态 关 系 存在 较 大 的 个 体 差 异 。 如 ， 对 于 个 体 模型 中 吸烟 欲望 的 自 回 
个 体 在 前 一 个 时 间 点 较 高 的 吸烟 欲望 均 会 预测 下 一 个 时 间 点 更 高 的 吸 类 
较 多 个 体 模型 中 二 者 不 存在 预测 作用 (47%)， 甚 至 可 能 存在 负 向 的 自 回 
在 9 个 个 体 模型 中 ,吸烟 欲望 可 以 预测 下 一 时 间 点 的 抑郁 水 平 ， 其 中 4 人 为 负 向 预测 ，5 人 


为 正 向 预测 。 


因 欲 望 和 抑郁 水 平 之 间 的 动 


归 效应 Purge Surgelag’ 较 多 
因 欲 望 (47%)， 但 也 有 
归 路 径 (6%)。 另 外 ， 
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5 个 体 模型 (随机 示例 ) 


YE: urge: 吸烟 欲望 ，dep: 抑郁 水 平 。 红 线 表 示 正 向 的 效应 ， 蓝 线 表 示 负 向 的 效应 ， 实 线 代 表 同 时 效应 ， 虚 线 表 示 清 后 


效应 (包括 自 


回归 和 交叉 滞后 效应 )。 


4.3.4 小 结 


本 节 分 别 利 用 MLM, RDSEM 和 GIMME 进行 实证 演示 和 分 析 。RDSEM 以 DSEM 作 


为 基础 ， 


因此 本 节 也 较 好 展示 了 DSEM 的 分 析 过 程 及 思路 。MLM 和 DSEM 均 为 自 上 而 下 


的 建 模 方 法 ， 相 比 MLM, DSEM 在 分 析 密 集 追 踪 数 据 时 存在 诸多 优势 : (1) 可 以 方便 地 对 变 


量 间 的 自 


个 因 变 量 


回归 和 交叉 滞后 效应 进行 分 析 ( 包 括 相 反方 向 的 交叉 滞后 效应 ); (2) 可 以 同时 考虑 多 


; (3) 对 模型 随机 效应 的 数量 几乎 没有 限制 ，(4) 可 以 在 水 平 1 输出 不 同 个 体 中 各 项 
标准 化 系数 的 平均 结果 ; (5) 采 用 潜 变 量 中 心 化 ，(6) 可 以 考虑 残 差 项 方差 的 随机 效应 。 如 果 
研究 者 希望 更 加 系统 、 深 入 地 学 习 D 


SEM 在 Mplus 中 的 代码 实现 ， 可 以 参考 McNeish 和 


Hamaker(2019) 的 文章 。 
与 MLM 和 DSEM 不 同 , GIMME 是 一 种 典型 的 自 下 而 上 的 建 模 方法 。 与 DSEM 相 比 ， 
GIMME 的 优势 主要 体现 在 它 可 以 单独 输出 各 个 个 体 模型 的 结果 ， 人 允许 不 同 个 体 模型 的 结构 


存在 异 质 


性 ， 从 而 更 有 针对 性 地 呈现 每 个 个 体 中 变量 相互 作用 的 模式 。 然 而 ，GIMME 不 能 


将 不 随时 间 变 化 的 变量 纳入 模型 ， 因 此 不 能 方便 地 考虑 个 体 特异 的 变量 (如 工作 压力 、 家 庭 


压力 ) 对 个 体内 部 变量 间作 用 关系 的 影响 。 


4.4 其 他 模型 


除了 以 上 介绍 的 模型 方法 ,在 密集 追踪 数据 中 也 存在 其 他 模型 方法 的 应 用 。 如 ,考虑 到 
心理 学 研究 中 常常 出 现 潜 变量 的 情况 ， 研 究 者 们 将 MLM 和 SEM 相 结合 ， 提 出 多 层 结构 方 
程 模型 (Multilevel Structural Equational Model, MSEM; Gottfredson et al., 2009; Johnson et al., 
2005). 5 MLM 类 似 ，MSEM 在 探究 个 体内 部 的 动态 过 程 以 及 个 体 特异 性 等 方面 也 存在 较 
大 的 局 限 性 。 如 研究 者 希望 更 深入 地 了 解 MSEM, 可 以 参阅 Rabe-Hesketh 等 (2007), Preacher 


等 (2010) 以 及 Holtmann 等 (2016) 的 文章 。 
多 层 VAR 模型 (Multilevel Vector Autoregressive Model, mlVAR; Bringmann et al., 2015) 则 
是 在 VAR 模型 的 基础 上 ， 进 一 步 允 许 对 多 个 个 体 的 密集 追踪 数据 进行 分 析 ， 也 是 目前 较为 


常见 的 用 于 密集 追踪 的 方法 之 一 (Bringemann et al., 2016; Jongerling et al., 2015; Piccirillo & 


Rodebaugh, 2019)。 特别 是 其 与 网 络 分 析 相 结合 ， 以 及 其 R 包 mIVAR 的 提出 (Bringmann etal., 


j 2013; Epskamp et al., 2016), #74 mIVAR 的 应 用 相对 更 加 简便 ， 且 同样 可 以 输出 与 GIMME 
= 类 似 的 模型 网 络 ， 便 于 在 临床 上 识别 症状 网 络 中 的 症状 中 心 (Symptom Centrality; Bringmann 
etal., 2015)。 然 而 ,与 GIMME 不 同 ，mlVAR 本 质 上 采用 的 是 基于 整体 的 建 模 思想 。mlVAR 

有 具有 测量 时 间 点 的 数量 要 求 较 低 、 适 用 于 高 维 数据 分 析 的 优点 , 但 它 只 能 进行 显 变量 的 建 模 ， 


且 难 以 对 复杂 模型 进行 估计 (Bringmann et al., 2015; Epskamp et al., 2018)。 由 于 篇 幅 有 限 ， 本 


© 文 仅 对 mlVAR 进行 简单 介绍 ， 更 详细 的 介绍 可 参见 Bringmann 等 人 (2015) 的 文章 。 
另外 ， 考 虑 到 随 着 时 间 的 推移 ， 变 量 动 态 变化 过 程 中 可 能 存在 一 定 的 时 变 效 应 (Time- 


Varying Effect), Bringmann 等 人 (2018) 在 VAR 的 基础 上 进一步 将 时 变 效 应 纳入 考虑 , 提出 时 


| 


变 向 量 自 回 归 模 型 (Time-Varying Vector Autoregression, TV-VAR; Bringmann et al., 2018). TV- 
VAR 最 初 起 源 于 计量 经 济 学 领域 , 但 近年 来 也 被 应 用 于 分 析 心 理 过 程 的 动态 变化 (Bringmann 
et al., 2017; Bringmann et al., 2018; Chow et al., 2017; Haslbeck & Waldorp, 2020)。 与 前 面 介绍 
的 模型 方法 不 同 , TV-VAR 不 需要 假设 时 间 序 列 的 稳定 性 (Stationarity), 而 是 允许 模型 参数 随 
时 间 发 生变 化 ( 唐 文 清 等 , 2020; Bringmann et al., 2018)。 然 而 ，TV-VAR 对 于 观测 数据 点 的 数 
量 要 求 较 高 (一 般 需 要 对 每 个 被 试 测量 超过 100 个 时 间 点 )， 且 目前 只 能 用 于 两 变量 的 建 模 分 


析 (Piccirillo & Rodebaugh, 2019)。 感 兴趣 的 研究 者 可 以 参阅 唐 文清 等 人 (2020) 以 及 Bringmann 


等 人 (2018) 的 文章 。 


4.5 密集 追踪 数据 分 析 方 法 的 选择 策略 


为 了 帮助 研究 者 针对 自身 研究 的 实际 情况 , 灵活 选择 适合 的 模型 方法 , 这 里 对 各 种 主流 


模型 方法 的 优 缺 点 及 选择 策略 进行 归 类 和 总 结 ， 以 供 参考 (Asparouhov etal., 2018; Gates etal., 


2020; Gates & Molenaar, 2012; Piccirillo & Rodebaugh, 2019): 


表 1 密集 追踪 的 主流 方法 优 缺 点 


模型 方法 优点 缺点 
假设 不 同 个 体 间 变量 的 相互 作用 机 制 是 同 质 
MLM/MSEM ”简单 易 行 ， 对 时 间 点 数量 要 求 相对 较 低 。 的 ; 难以 分 析 个 体内 部 的 动态 变化 过 程 及 机 


对 时 间 点 数量 要 求 相对 较 低 ; 适 


mlVAR 


适用 于 较 多 


变量 间 的 动态 交互 ， 结 果 可 视 性 强 


可 以 将 随机 效应 分 解 为 个 体 


和 时 间 两 个 


方面 的 来 源 ; 在 时 间 节 点 和 随机 效应 的 数 


DSEM 


IN 
EA] 

mal 
可 


量 上 限制 较 少 ; 允许 测量 时 间 节 点 之 间 间 


制 ， 建 模 不 够 灵活 


假设 不 同 个 体 间 变量 的 相互 作用 机 制 是 同 质 
的 ; 不 考虑 潜 变 量 因 子 结构 


假设 不 同 个 体 间 变量 的 相互 作用 机 制 是 同 质 
的 ， 需 要 对 被 试 测量 较 多 时 间 点 


针对 单个 被 试 进行 建 模 , 考虑 个 体内 部 的 


(LV-)GIMME ”动态 变化 过 程 及 机 制 的 异 质 怕 


Es 适用 于 较 


多 变量 间 的 动态 交互 ， 自 动 化 搜索 


可 能 出 现 模 型 过 拟 合 ， 无 法 将 几乎 不 随时 间 
变化 的 变量 纳入 模型 ， 需 要 对 单个 被 试 测量 
较 多 时 间 点 


a 


YE: GIMME: 组 迭代 多 模型 估计 (Group Iterative Multiple Model Estimation); LV-GIMME: 基于 洪 变 量 进 行 建 模 的 


GIMME(Latent Variable GIMME); MLM: 多 层 线性 模 


型 (Multilevel Model); MSEM: 多 层 结构 方程 模型 (Multilevel Structural 


Equational Model); mlVAR: 多 层 VAR 模型 (Multilevel VAR); DSEM: 动态 结构 方程 模型 (Dynamic Structural Equation 


Model). FE]. 
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5 讨论 


密集 追踪 数据 分 析 


6 密集 追踪 数据 方法 的 选择 策略 流程 图 


密集 追踪 方法 对 个 体 展开 多 个 密集 时 间 节 点 的 数据 奶 踪 ,能 更 好 地 研究 个 体 的 变化 过 程 


和 内 部 机 


出 ， 在 外 部 效 度 、 数 据 分 析 方 法 等 方面 相 比 传统 的 追踪 研究 有 具有 较 大 优势 。 为 了 鼓 


励 更 多 的 研究 者 在 密集 追踪 数据 中 关注 动态 变化 的 过 程 、 不 同时 间 节 点 间 存 在 的 自 回 归 效 应 


等 ， 本 文 主要 介绍 了 两 种 密集 追踪 模型 方法 的 新 


运用 和 


模 思 想 


Cun 


考虑 (Asparouhov et al., 2018)， 而 GIMME 从 纯粹 的 基于 个 人 的 分 析 方 法 
单个 个 体 的 模式 ; 
异性 研究 的 1 


展 ， 和 希望 能 够 增进 研究 者 们 对 其 的 理解 、 
步 发 展 。 相 比 传统 的 分 析 方 法 , DSEM 和 GIMME 都 更 多 地 纳入 了 时 间 序 列 的 建 


并 更 加 强调 对 于 个 体 间 差 异 的 关注 。DSEM 允许 研究 者 对 任意 数量 的 随机 效应 进行 


8 发， 更 直接 地 对 
行 探索 。 这 些 分 析 方 法 对 个 体 差 异 的 高 度 关注 , 反映 了 近年 来 对 于 个 体 特 


昌 导 和 发 展 的 趋势 。 


值得 注意 的 是 ，GIMME 的 一 大 特点 在 于 其 模型 的 建立 是 探索 性 或 者 半 验 证 性 的 ， 这 与 


传统 心理 学 


究 中 基于 前 人 研究 及 已 有 的 知识 经 验 提出 假设 的 模型 、 利 用 大 规模 的 数据 收集 


对 模型 i 


动 的 方式 对 变量 间 的 关系 进行 自由 探索 , 可 能 可 以 避免 先入 为 主 的 固化 思维 ,为 研究 者 提供 


行进 


步 验 订 


E 的 思想 存在 很 大 的 不 同 。 这 种 探索 性 的 建 模 思路 人 允 询 


F 研 究 者 以 数据 驱 


更 多 新 的 发 现 和 启发 。 然而， 探索 性 分 析 常 常 也 伴随 模型 过 拟 合 、 结 果 难 以 解释 等 问题 ， 未 
来 还 需要 研究 者 进一步 探究 对 其 停止 搜索 规则 的 优化 设置 。DSEM 对 模型 中 随机 效应 的 数 


量 几 乎 没有 限制 ,这 在 改善 模型 拟 合 的 同时 也 可 能 导致 模型 过 于 复杂 、 无 法 泛 化 到 更 广泛 群 
体 的 问题 , 为 了 尽 可 能 减少 可 重复 性 方面 的 问题 , 建议 研究 者 在 实际 研究 前 对 其 实验 设计 及 
分 析 计划 进行 预 注 册 ， 并 在 文章 中 尽 可 能 详尽 地 报告 其 数据 收集 、 处 理 及 分 析 等 过 程 。 

另外 , 目前 对 于 密集 追踪 数据 的 动态 分 析 大 多 假设 测量 时 间 点 是 离散 化 的 , 这 种 假设 在 
连续 时 间 测 量 ， 以 及 测量 时 间 间 隔 高 度 随 机 的 情况 下 可 能 存在 一 定 的 问题 (Hamaker et al., 
2018)。 连 续 时 间 建 模 (Continuous Time Modeling) 可 以 较 好 地 解决 时 间 离 散 化 假设 带 来 的 问 
题 ， 但 它 需 要 涉及 复杂 的 微分 计算 (Voelkle et al., 2018)， 目 前 应 用 研究 者 对 于 如 何在 实证 分 
析 中 实现 连续 时 间 建 模 还 不 够 熟悉 。 未 来 还 需要 更 多 研究 者 进一步 探索 连续 时 间 模 型 的 简化 
计算 与 自动 化 实现 ， 为 应 用 研究 者 提供 更 灵活 的 建 模 选 择 。 

密集 追踪 研究 关注 个 体内 行为 、 认 知 、 心 理 等 的 动态 变化 过 程 ， 这 些 内 在 过 程 及 机 制 往 
往 十 分 复杂 ， 可 能 存在 非 线性 的 作用 模式 (如 ， 两 变量 间 存 在 交互 作用 )。 目 前 密集 追踪 数据 
较为 主流 的 分 析 方 法 大 多 仅 考虑 线性 效应 的 情况 (如 ，DSEM 和 GIMME)， 但 仅 通 过 线性 关 


系 对 这 些 动态 效应 进行 解释 很 可 能 是 不 够 的 。 最 近 , Kelava 和 Brandt(2019) 将 DSEM 与 动态 


lin 


潜 类 别 分 析 (Dynamic Latent Class Analysis; Asparouhov et al., 2017) 相 结合 ， 提 出 非 线性 的 动 

态 潜 类 别 结构 方程 模型 (Nonlinear Dynamic Latent Class Structural Equation Modeling, NDLC- 
SEM)， 人 允许 在 组 间 和 组 内 水 平 进行 非 线 性 效应 的 分 析 (Kelava & Brandt, 2019)。 然 而 ， 目 前 
允许 对 密集 追踪 数据 进行 非 线 性 建 模 的 模型 方法 仍 较为 有 限 , 其 在 不 同情 境 下 的 具体 表现 还 
不 明晰 ， 未 来 需要 更 多 方法 学 研究 者 对 其 进行 更 进一步 的 研究 。 

密集 追踪 研究 的 数据 收集 往往 基于 被 试 的 真实 生活 场景 ， 具 有 较 高 的 外 部 效 度 。 然 而 ， 
真实 的 测量 环境 往往 也 更 加 复杂 , 研究 者 可 能 需要 相应 增加 对 潜在 干扰 因素 的 动态 测量 , 并 
将 其 纳入 模型 进行 分 析 。 此 时 ， 模 型 中 涉及 的 变量 数目 较 多 ， 因 此 采用 GIMME 对 数据 进行 
分 析 可 能 是 更 合适 的 。 然 而 ，GIMME 对 于 单个 个 体 的 测量 时 间 点 数 要 求 较 高 ， 这 大 大 增加 
了 应 用 研究 者 密集 追踪 数据 的 收集 难度 。 未 来 可 以 进一步 探索 适用 于 高 维 密集 追踪 数据 、 且 
对 测量 时 间 点 要 求 相对 较 低 的 模型 方法 , 以 及 进一步 发 展 密集 追踪 数据 的 收集 方法 , 尽 可 能 
降低 其 收集 的 难度 。 
总 之 , 密集 追踪 方法 的 一 大 优势 在 于 探究 个 体内 部 的 动态 变化 过 程 及 其 机 制 。 虽然 目前 
密集 追踪 数据 分 析 在 方法 学 领域 已 有 较 多 新 的 进展 , 但 它 在 实际 研究 中 应 用 不 足 , 许多 应 用 
研究 者 缺乏 对 于 时 间 这 一 要 素 的 关注 (Cronin & Vancouver, 2019; Voelkle et al., 2018; Zhou et 
al.,2020)。 现 在 ， 较 多 密集 追踪 数据 的 分 析 方 法 已 经 可 以 方便 地 通过 软件 实现 (如 ，DSEM 和 
GIMME)， 大 大 减少 了 实证 研究 应 用 的 阻碍 。 和 希望 本 文 可 以 帮助 实证 研究 者 们 了 解 、 掌 握 更 


+ 


4/ 


E 


es 


灵活 的 建 模 和 分 析 方 法 ， 更 好 地 发 挥 密集 追踪 数据 的 优势 ， 尽 可 能 减少 对 数据 结果 的 误 读 。 
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Intensive longitudinal data analysis: Models and application 
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Abstract: In the fields of psychology, education, and clinical science, researchers have devoted 
increased attention to the dynamic changes and personalized modeling of individuals’ behaviors, 
minds, and treatment effects over time. Intensive longitudinal data is a set of measures collected at 
multiple time points with higher frequency over shorter periods. Thus, it can be used in the analysis 
of the dynamics and mechanisms of within-person processes. In recent years, intensive longitudinal 
design has become one of the most prominent and promising approaches in psychological research. 
However, many of these researches still rely on traditional data analysis methods. Many models 
have been proposed to analyze intensive longitudinal data, including top-down approaches (e.g., 
dynamic structural equation model, DSEM) and bottom-up approaches (e.g., group iterative 
multiple model estimation, GIMME). Both of the methods can conveniently model autoregressive 
and cross-lagged effects in intensive longitudinal data. 
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附录 
附录 A: MLM 与 RDSEM 实证 演示 Mplus 代码 


1、MLM 代码 : 
DATA: 
FILE = [DATA LOCATION PATH)]\ two-level data with Trend.csv; 
VARIABLE: 
NAMES ARE URGE DEP JS HS PERSON TIME; 
USEVARIABLES URGE DEP JS HS TIME; 
cluster=PERSON; 
BETWEEN = JS HS; 
WITHIN = TIME DEP;! Time is only used in the Within-Level model; 


DEFINE: 
CENTER DEP (groupmean); 


ANALYSIS: 


TYPE = TWOLEVEL RANDOM;! Two-Level Model with paths that are latent variables; 


MODEL: 
“%WITHIN% 
bli| URGE ON DEP; 
! urge to smoke is regressed on Lag-1 urge to smoke 


trend | URGE ON TIME:! urge is regressed on time to capture the trend over time, the !slope 


is latent; 


%BETWEEN% 
URGE bli; 
URGE with bli TREND;bli with TREND; 


[URGE]; !mean intercept of urge to smoke; 
[bli]; !mean of the regressive slope; 
[TREND]; 


URGE; !urge to smoke intercept variance; 
bli; !regressive slope variance; 
TREND; 


URGE on JS HS; 

!The person-specific intercept is predicted by Job Stress and Home Stress; 
blion JS HS; 

!The regressive slope is predicted by Job Stress and Home Stress; 
TREND on JS HS; 


OUTPUT: TECH1I TECH8; 


2. RDSEM 代码 : 

DATA: 
FILE = [DATA LOCATION PATH]\ two-level data with Trend.csv; 

VARIABLE: 
NAMES ARE URGE DEP JS HS PERSON TIME; 
USEVARIABLES URGE DEP JS HS TIME; 
Lagged = URGE(1) DEP(1);!Created a Lag-1 variable for urge to smoke and depression; 
CLUSTER=PERSON; 
BETWEEN = JS HS;! Job Stress and Home Stress are only used in the Between-Level model; 
WITHIN = TIME;! Time is only used in the Within-Level model; 


ANALYSIS: 
TYPE = TWOLEVEL RANDOM; ! Two-Level Model with paths that are latent variables; 
ESTIMATOR = BAYES;!Lagged variables can only be estimated with Bayes in Mplus; 
BITERATIONS = (1000);!Run at least 1000 iteration of the MCMC algorithm; 
BSEED=1028; 


MODEL: 
%WITHIN% 
PHI1 | URGE’ ON URGE“1; 
! urge to smoke is regressed on Lag-1 urge to smoke 
PHI? | DEP’ ON DEP^I; 
! depression is regressed on Lag-1 depression 
PHI3 | DEP’ ON URGE“ 1; 
! Depression is regressed on Lag-1 urge to smoke 
PHI4 | URGE’ ON DEP“1; 
! urge to smoke is regress on Lag-1 depression 
TREND! | URGE ON TIME:! urge is regressed on time to capture the trend over time 
TREND2 | DEP ON TIME;! dep is regressed on time to capture the trend over time 
LogV1 | Urge; 
! the Within-Level residual variance of urge to smoke 
LogV2 | DEP; 
! the Within-Level residual variance of depression 


%BETWEEN% 

[URGE]; !mean intercept of urge to smoke 
[DEP]; !mean intercept of depression 

[PHI1]; !mean of urge autoregressive slope 
[PHI2]; !mean of depression autoregressive slope 
[PHI3]; !mean of dep on urge&1 slope 

[PHI4]; !mean of urge on dep&1 slope 


[TREND 1]; 

[TREND2]; 

[LogV 1]; 

llog of the mean of the urge residual variance 
[LogV2]; 

!log of the mean of the urge residual variance 


URGE; !urge to smoke intercept variance 

DEP; !depression intercept variance 

PHI1; !urge to smoke autoregressive slope variance 

PHI2; !depression autoregressive slope variance 

PHI3; !dep on urge&1 slope variance 

PHI4; !urge on dep&1 slope 

TREND]; 

TREND2， 

LogV1; 

!between-person variance of urge to smoke residual variance 


LogV2; !between-person variance of depression residual variance 


URGE on JS HS; 

DEP on JS HS; 

!The person-specific intercepts are predicted by Job Stress and Home Stress; 
PHIl1 on JS HS; 

PHI2 on JS HS; 

PHI3 on JS HS; 

PHI4 on JS HS; 

!The person-specific slopes are predicted by Job Stress and Home Stress; 
LOGV1 on JS HS; 

LOGV2 on JS HS; 

!The person-specific residual variances are predicted by Job Stress and Home Stress; 


URGE DEP PHI1-PHI4 LOGV1 LOGV2 WITH URGE DEP PHI1-PHI4 LOGV1 LOGV2; 


OUTPUT: TECH1 TECH8 STDYX; 


附录 B: GIMME 实证 演示 R 代码 


library(DiagrammeR) 
library(gimme) 
library(mvtnorm) 
library(dplyr) 


library(xlsx) 


#1、 数 据 预 处 理 
wd<-"~""# 设置 文件 夹 路 径 


Setwd(wd) 


dat<- read.csv("Two-Level Data with Trend.csv",header=F) # 读 取 数据 文件 


colnames(dat)<-c("URGE", "DEP", "JS", "HS", "PERSON", "TIME") # 指定 各 变量 名 ， 其 中 


URGE 指 吸烟 欲望 ，DEP 指 抑郁 水 平 


datawd<-paste0(wd,"/data") 

setwd(datawd) 

for(i in 1:length(unique(datSPERSON))) { 
# 选取 个 体 i 的 数据 进行 预 处 理 


person dat<-dat[which(dat$PERSON==i),c("URGE","DEP","TIME")] 


# 分 别 去 除 变 量 URGE 和 DEP 随时 间 变 化 的 整体 趋势 (detrend) 
temp = lIm(person dat$URGE~person dat$TIME) 
person_dat$urge<-residuals(temp) 

temp2 = lm(person_dat$DEP~person_dat$TIME) 


person_dat$dep<-residuals(temp2) 


# 将 个 体 i 的 数据 单独 写 入 excel 文档 "data i.csv" 


write.csv(person_dat[,c("urge","dep")], 


paste("data",i,".csv"), row.names = FALSE) 


#2. NJH gimmeSEM 函数 输出 群体 和 个 体 模型 结 

fit<-gimmeSEM( 
data =datawd, # 存放 数据 的 文件 夹 路 径 
out = pasteO(wd,"/output"), # 存放 输出 结果 的 文件 夹 路 径 
sep=".", # 声明 数据 分 隔 符 ， 这 里 cov 格式 的 数据 分 隔 符 为 逗号 


是 否 包括 变量 名 


oO 


header=T, # 声明 数据 


202009.00002v2 


chinaXiv 


ar=TRUE, # 指定 模型 中 是 否 包括 自 回 归 ， 默 认为 TRUE 

plot=TRUE, # 指定 结果 是 否 输出 模型 图 ， 默 认为 TRUE 

groupcutoff = .75 # 路 径 可 以 被 纳入 群体 模型 的 临界 值 (显著 改善 个 体 模型 的 最 小 
比例 ， 默 认为 75%) 
) 


